PARTIAL DIFFERENTIAL OPERATORS WITH NON-NEGATIVE 
CHARACTERISTIC FORM, MAXIMUM PRINCIPLES, AND UNIQUENESS 
FOR BOUNDARY VALUE AND OBSTACLE PROBLEMS 

PAUL M. N. FEEHAN 

Abstract. We prove weak and strong maximum principles, including a Hopf lemma, for classical 
solutions to equations defined by linear, second-order, partial differential operators with non- 
negative characteristic form (degenerate-elliptic operators), in the presence of a second-order 
boundary condition of Ventcel type along the degeneracy locus of the principle symbol of the 
operator on the domain boundary. We apply these maximum principles to obtain uniqueness 
and a priori maximum principle estimates for classical solutions to boundary value and obstacle 
problems defined by these degenerate-elliptic operators, again in the presence of a second-order 
boundary condition, for Dirichlet or Neumann boundary conditions along the complement of the 
degeneracy locus. We also prove weak maximum principles and uniqueness for solutions to the 
corresponding variational equations and inequalities defined with the aide of weighted Sobolev 
spaces. The domain is allowed to be unbounded when the operator coefficients and solutions 
obey certain growth conditions. 
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1. Introduction 

The classical maximum principles of Fichera |34l [35] and Olemik, and Radkevic |58l [63l [64] 
provide uniqueness theorems for degenerate elliptic and parabolic boundary value problems which 
do not take into account a more modern view of the appropriate function spaces in which unique- 
ness is sought, such as [iHl [HI [20l [26l [271 EH 133 [47]. Indeed, their maximum principles lead to 
the imposition of additional Dirichlet boundary conditions which are not necessarily motivated 
by the underlying application, whether in biology, finance, or physics. These additional Dirich- 
let boundary conditions, usually for certain ranges of parameters defining the partial differential 
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equation, are often less natural than the physically-motivated regularity properties suggested 
by choices of appropriate weighted Holder [19l [201 123 or Sobolev spaces [18l EH ST] , which 
automatically encode special regularity or integrability up to portions of the domain boundary 
where the operator becomes degenerate. In the case of weighted Holder spaces, these boundary 
regularity properties may be viewed as special cases of second-order or Ventcel boundary con- 
ditions [H El [68] (see Remark 12. 5p . However, the question of exactly how regular the solution 
should be near these boundary portions is delicate: asking for too much regularity (such as 
up to the boundary) and the boundary value problem may not be well-posed, while asking for too 
little regularity (such as up to the boundary) may lead to the requirement of an unphysical 
Dirichlet boundary to ensure the problem is well-posed, with the unintended consequence that 
the solutions thus selected can be no more than continuous up to the boundary; an illustration 
of this point is provided by Example 11.31 

1.1. Formulation of boundary value and obstacle problems for linear, second-order 
partial differential operators with non-negative characteristic form. Let ^ C M'', d>2, 
be a domain, that is, a connected, open, possibly unbounded subset with boundary dff and 
suppose that S ^ di^ is a relatively open subset. We shall consider the question of uniqueness 
and a priori estimates for the boundary value problem, 

Au= f on ^, 

(1.1) 

u = g on dff\T,, ^ ^ 

and for the obstacle problem, 

min|j4ti — f,u — th} = a.e. on ^, 

(1.2) 

u = g on dff\j:, ^ ^ 

where ip < g on dG \ E, and 

Au := - tv{aD^u) - {b, Du) + cu, ue C°°(^), (1.3) 
where a: G ^ W^'^'^ has the property thalQ a(x) is a non-negative definite matrix for all x G i^, 

(a(x)r/,r/) > 0, Vx G ^, V?] G M'^. 

and choos^ 

S = int ({x G diff : a{x) = 0}) , 

where int(S') denotes the interior of a subset S of a topological space. Throughout this article 
we shall allow S to be non-empty. Finally, b : G ^ is a vector field and c : — >■ M is a 
non-negative function. We use D^u and Du to denote the Hessian matrix and gradient of n, 
respectively, so tr(aL'^n) = a^^Uxix^ and {b,Du) = b^Uxi, where Einstein's summation convention 
is used throughout this article. Following [58], we call A in (II. 3p a linear, second-order partial 
differential operator with non-negative characteristic form. Such operators are more commonly 
known as "degenerate elliptic" when S is non-empty (see Remark 14.41 for definitions of elliptic, 
strictly elliptic, and uniformly elliptic operators). 



^Though unnecessary for the maximum principle, an additional assumption of symmetry ensures that we may 
find a square root, a : 6 ^ R'**^'' such that a = ^oa* and consider the martingale problem for (a, 6) |67) . 

For the sake of notational consistency in this article, we do not always distinguish between situations where 
E ^ d6 may be any open subset rather than the zero locus of the matrix a on dG and instead rely on the context 
to make the distinction clear. 
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For the sake of clarity, we shall confine our attention to linear operators in this article, although 
many of the results can be seen to extend to semilinear operators, 

S{u) = - tr{aD^u) - {b, Du) + c(-, n), n G C'^{ ff), 

or quasilinear operators [iOl §10], 

Q{u) = - ti{aD'^u) - b{-,u, Du), u £ C'^{ff), 

as well as linear and quasilinear parabolic operators with non-negative characteristic form. We 
shall discuss such situations in separate articles. 

The conditions on the coefficients of A in (11.3P assumed in this article are mild enough that 
they allow for many examples of partial differential operators. A, with non-negative characteristic 
form and which are of interest in mathematical biology, finance, and physics. Before proceeding 
to a discussion of our main results, we shall first provide some specific examples of operators to 
which our results apply. 

1.2. Examples. We begin by describing a family of examples which includes certain stochas- 
tic volatility models occurring in mathematical finance [l5] and the linearization of the porous 
medium equation [19] . 

Example 1.1 (Affine coefficients and degeneracy on the boundary of a half-space). Suppose the 
coefficients of A in (jl.Sp are affine functions of x E M'^, with a[x) positive definite for all x E H, 
where HI = M.^~^ x M+ is a half-space and M+ = (0, oo), while a{x) = if x E dM. Then 

Au = -trixaaiD^u) - {bo + bix,Du) + {co + {ci,x))u, ueC°°{M), (1.4) 

where oi, 5i E M'^^'^ and bo, c\ E M*^ and cq E M. Thus, A is an elliptic partial differential operator 
on C°°(H) which becomes degenerate along the boundary S = 5BI = \xd = 0} of the half-space 
M = {xrf > 0}. 

When a\ is symmetric, the operator —A is the generator of a degenerate diffusion process. 
Imposing the condition (6, n) > along 5IHI, where n is the inward-pointing unit normal vector 
field, ensures that the diffusion process remains in the half-space H = {x^ > 0} if started in H. 
Since n = e^, this translates to the requirement that b'^{x) > for all x E dM, and thus 6q > 
and bf = for 1 < j < d - 1 and bf > 0. 

To ensure uniqueness of solutions to (II. ID on H via our maximum principle, it is necessary 
(though not sufficient) to impose the condition c > on H, and hence cq > and c\ = for 
1 < i < d - 1 and cf > 0. 

Examples of this kind occur frequently in mathematical finance, such as the Heston stochastic 
volatility process [45], where a{x) = x^ai for x E H, and oi E M'^^'^ is positive definite, and d = 2. 
See Example 11.21 for a description of the Heston process generator and an important example of a 
degenerate affine process. The linearization of the porous medium operator is another important 
example of this type; see Example 11.41 for details. □ 

Example 1.2 (Elliptic Heston operator). The generator of the Heston process |45j provides a 
well-known example in mathematical finance of the operator (jl.4p when d = 2: 

Au := - Y {uxix2 + '^Q(7Ux^x2 + o-^^ixaxa) - {r - q - X2/2)ux^ - k{9 - X2)u^2 + '^^^ (1-5) 

where q > 0,r > 0,k > 0,^ > 0,0" > 0, and q ^ ( — 1,1) are constants (their financial interpreta- 
tion is provided in [l5]), and u E C°°(]H), with EI = M x M+. □ 
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Example 1.3 (Elliptic Feller square-root or Cox-Ingersoll-Ross operator). The generator of the 
Feller square-root process [33], known as the Cox-Ingersoll-Ross process in mathematical finance 
[13], [11], [Ml Example 6.5.2], provides a simple example of the operator (|1.4p when d = 1: 

Au := -—xuxx - f^iO - x)ux + ru, (1.6) 

where u € C°°(M4.), and which takes the form, after a change of variables, of the Kummer 
operator 

Bv := -xvxx -{13- x)vx + av, (1.7) 
where v G C°°(M+) and /3 := 2k^/(T^ > and a := r/n > 0. The homogeneous Kummer 
equation, Bv = on M+, has two independent solutions, the confluent hypergeometric function 
of the first kind M{a, x) (or Kummer function) and the confluent hypergeometric function of 
the second kind U{a,l3;x) (or Tricomi function) [21 §13.1.2 & §13.1.3]. The solution M(a, /3;x) is 
in Ci^,(t+), with M(a, /3; 0) = 1, A-4(q,/3;0) = a//3, and A4^(a,/3;0) = a(a + l)/(/3(/3 + 1)) [21 
§13.4.9]. Near x = 0, the solution f7(a, /3;x) is comparable to x^~^ when /3 7^ 1 and log a; when 
^ = 1 [21 §13.5.6-12], and thus is in Ci°c(^+) < /3 < 1, with Uxia,(3;x) comparable to x"'^ 
and Uxx{ci,P',x) comparable to x~^~^ near x = 0, for any /3 > [21 §13.4.22], and thus not even 
in C/„^(t+) when /3 > 0. 

Example 1.4 (Linearization of the porous medium operator). In a landmark article, Daskalopou- 
los and Hamilton [19] proved existence and uniqueness of C°° solutions, u, to the Cauchy problem 
for the porous medium equation [191 P- 899] (when d = 2), 

d 

-ut + Y,i^"')^^x, = (0' X 0) = 5 on (1.8) 

i=l 

with constant m > 1 and initial data, g > 0, compactly supported on W^, together with C°°- 
regularity of its free boundary, d{u > 0}, provided the initial pressure function is non-degenerate 
(that is, Du"^~^ > o > 0) on boundary of its support at t = 0. Their analysis is based on their 
development of existence, uniqueness, and regularity results for the linearization of the porous 
medium equation near the free boundary and, in particular, their model linear degenerate operator 
\19\ p. 901] (generalized from d = 2 in their article), 

d 

Au := -XdY,^^,^, - nGC~(M), (1.9) 

i=l 

where /3 is a positive constant, analogous to the combination of parameters 2K9/a'^ in (jl.5p . and 
H = ]R.'^~^ X M_|_, following a suitable change of coordinates [19\ p. 941]. The same model linear 
degenerate operator (for d > 2), was studied independently by Koch [171 Equation (4.43)] and, 
in a remarkable Habilitation thesis, he obtained existence, uniqueness, and regularity results for 
solutions to (jl.Sp which complement those of Daskalopoulos and Hamilton [19] . □ 

Example 11.11 describes a class of elliptic differential operators which become degenerate along 
the boundary of a half-space, M.'^~^ x M+; their coefficients are affine functions of x G M"^ and 
M'^"^ X is a state space for the corresponding Markov process when the coefficient a{x) 
is symmetric. More generally, mathematical finance and biology provide examples of elliptic 
differential operators which become degenerate along the boundary of a "quadrant", W^~"^ x M™, 
and M'^"™ x M™ is a state space for the corresponding Markov process. Examples primarily 
motivated by mathematical finance include affine processes [H [I2l [I7l [22l [23l [24l [36l [37] , which 
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may be viewed as extensions of geometric Brownian motion (see, for example, [66])! the Heston 

stochastic volatihty process [l5], and the Wishart process [12l [38l HU |l2l |33] . Examples of this 
kind which arise in mathematical biology include the multi-dimensional Kimura diffusions and 
their local model processes [271 Equations (1.5) & (1.20)]. Another example along these lines is 
provided by the articles of Athreya, Barlow, Bass, Lavrentiev, and Perkins [UJ El [7] concerning 
generators of super-Markov chains. 

Maximum principles addressing these more complex situations where A has a non-negative 
characteristic form which becomes degenerate along a stratified space T, ^ dff, and in the broader 
sense that 



rather than S = int({a; G dff : a{x) = 0}, are discussed in our companion article |29j . 

1.3. Summary of main results and outline of our article. We shall leave detailed state- 
ments of our main results to the body of our article and simply provide a short outline of our 
article here to facilitate the reader seeking a particular conclusion of interest. Part [1] of our ar- 
ticle ^ ^ and ^ develops weak and strong maximum principles for operators on smooth 
functions and applications to boundary value and obstacle problems, while Part [2] of our article 
(^ and ^ develops weak maximum principles for bilinear maps and operators on functions in 
Sobolev spaces and applications to variational equations and inequalities. 

1.3.1. Weak and strong maximum principles for operators on smooth functions and applications 
to boundary value and obstacle problems. We say that u E C^(^) H C"^(^ U S) n Cioc(^) obeys a 
second-order boundary condition along T, defined by a principal symbol a G Cioc(i^U S;M'^^'^) if 



and A in (jl.3p has the weak maximum principle property on 0" UT, for this class of functions if 



The innovation here, of course, lies in the nature of the boundary condition for u along S (see Re- 
mark [2J3] for a discussion of equivalent first-order boundary conditions). When 0' is unbounded, 
the conditions on u are augmented by an assumption that sup^tt < oo. 

We shall find it very convenient to cleanly separate a discussion of conditions which may lead to 
A in ()1.3p having the basic weak maximum principle property, which are provided in ^from their 
applications for boundary value and obstacle problems which flow from this abstract property 
and which are discussed in sections [2l [3l and [H 

In ^ we consider applications of the weak maximum principle property to boundary value 
problems, including a comparison principle for subsolutions and super solutions and uniqueness 
for solutions to the Dirichlet boundary problem (Proposition 12.161) . a priori estimates for 
subsolutions, supersolutions, and solutions (Proposition 12.191) . ^^^^ an extension for the case of 
subsolutions which obey a growth condition on unbounded domains (Theorem 12. 20p . 

Section l3l contains application of the weak maximum principle property to obstacle problems. 
We develop a theory of continuous supersolutions to the obstacle problem, modeled on the corre- 
sponding theory for the Dirichlet boundary value problem in [401 §2.8 &: §6.3] for linear, second- 
order, uniformly elliptic partial differential equations on bounded domains. Our development has 
some independent interest since it is a theory of classical rather than viscosity solutions, unlike 
that of jl6j . We then develop applications of the weak maximum principle property to C^'^ so- 
lutions to obstacle problems, including a comparison principle for continuous supersolutions and 



S = {x e 5^ : deta(x) = 0} 



tr {aD'^u) G Cioc(^ U S) and tr {aD^u) =0 on S 
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uniqueness for solutions to the obstacle problem (Theorem I3.15P and a weak maximum princi- 
ple and a priori estimates for solutions and continuous supersolutions to the obstacle problem 
(Proposition EHI). 

In ^ we prove a strong maximum principle and develop its applications to boundary value 
problems with Neumann boundary conditions. We first prove a Hopf boundary point lemma (see 
Lemma l4.ip for operators A in ()1.3p which may become degenerate on 9^ using a novel choice 
of barrier functions. We then apply our version of the Hopf lemma to prove a strong maximum 
principle suitable for operators with non-negative characteristic form (Theorem 14. 6p and corre- 
sponding uniqueness results for solutions to equations with Neumann boundary conditions along 
9^ \ S (Theorem IMl and Corollary SH]) . 

Finally, in ^ we establish specific conditions on the coefficients (a, b, c) which ensure that the 
operator A in (II. 3p has the weak maximum principle property on S, initially for bounded 
functions on bounded domains (Theorem 15. ip . and then for bounded functions on unbounded 
domains (Theorem 15. Sp . 

1.3.2. Weak maximum principles for bilinear maps and operators on functions in Sobolev spaces 
and applications to variational equations and inequalities. We next consider variational equations 
and inequalities defined by bilinear maps, o, on weighted Sobolev spaces, H^{0',tv), defined by 
weight, to G C{ff)nL^{&) with ro > on i^, and a degeneracy coefficient, d e C\oc{^) with i? > 
on 0", where a Borel measurable function, u on 0", lies in H^{&,to) if 



Given a non-empty open subset S ^ di? (which may be arbitrary for now), we say that a obeys 
the weak maximum principle property on UT, for functions u € H^{0', It)) if 



where i/(J(^ U S, w) is the closure of C^(^U S) in H^{ff,xo) and C^(^U S) C C7°°(^) is the 
linear subspace of smooth functions which have compact support in i^U S, and u < on 80' \ S 
in the sense of H^{0, to) if n"*" G Hq{0 U S, to), where u'^ = max{u, 0}. When is unbounded, 
the conditions on u are augmented by an assumption that ess sup^u < oo. 

In ^ we consider applications of the weak maximum principle property, including the com- 
parison principle (Proposition 16. 7p and a priori estimates (Proposition 16.10"]) for H^{0,to) su- 
persolutions and solutions to variational equations. We also show that when a bilinear map a 
on H^{0,to) has a weak maximum principle property for subsolutions (on unbounded domains) 
which are bounded above, the property may extend to subsolutions which instead obey a growth 



condition (Theorem 16. 16p . 



Section [7] contains applications of the weak maximum principle property to variational inequal- 
ities. We prove uniqueness for solutions to variational inequalities defined by bilinear maps a on 
H^{0',to), a comparison principle for supersolutions and uniqueness for solutions to variational 
inequalities (Theorem I7.2p . and a priori estimates (Proposition 17. 9^ for H^{0',to) supersolutions 
and solutions to variational inequalities. 

Finally, in ^ we prove the weak maximum principle for functions in weighted Sobolev spaces. 
We first show that a bilinear map of the form. 




a{u,v) < 0,Vt; G H^{0UT.,to) with v>0 a.e. on 0, 
u < on d0\T, in the sense of H^{0, to) 





a{u,v) = / (a* 
7^ 



UxiVxj + d^uvxj - h'uxiV + cuv) to dx, Vii, v G C^i^) 
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has the weak maximum principle property on ^ U S (Theorem 18. lip when G is bounded and its 
coefficients obey, for some positive constant, a.e. on G and all 77 G M'^, 

^\r]^ < {a7],r]) < Kt?!?]^, 

{b,v)<mvl 
\c\ < Kd, 

together with a non-negativity condition, 

a{l,v) > 0, Vt; G Hl{GU'L,'m) with u > a.e. on 

and a key weighted Sobolev inequality holds: given constants 1 < p < q < 00, there is a positive 
constant C such that, for any u G Li{G,tt)) with Du G LP{0','dto;R'^), one has 

Given additional constraints on the coefficients, we then extend our weak maximum principle to 
the case of bounded H^{G,m) functions on unbounded domains (Theorem IS.lSp . 

Our previous definition of the degeneracy locus, S, for a{x) is not meaningful when the function 
a: ^ M'^^'^ is only Borel measurable but, with the aid of the degeneracy coefficient, i?, it may 
be reinterpreted as 

S = int ({x G dG : i9(x) = 0}) , 

and, we then obtain the required weighted Sobolev inequality (Corollary 18. 6p for specific choices 
of t?(x), including a power of the distance, dist(x, S), from x G 1^ to the degeneracy locus, E ^ dG 
via a result of Koch [47] . 

Several appendices contain additional results and comments which shed further light on our 
conclusions. Appendix [A] provides examples illustrating when the weak maximum principle holds 
for functions obeying suitable growth conditions on unbounded domains. Appendix iBl contrasts 
the weighted Sobolev inequality (Corollary 18. 6p due to Koch with one due to V. Maz'ya- In 
appendix [Cl we compare the weak maximum principles and uniqueness theorems provided by 
our article with those of Fichera, Olemik, and Radkevic [63] in the case of the elliptic Heston 
operator, in Example 11.21 on a domain G and show that those of Fichera, Olemik, and 
Radkevic are strictly weaker. 

1.4. Comparison with previous research. A detailed comparison between our weak maxi- 
mum principle for classical subsolutions (Theorems 15. II and l5.3p and that of Fichera |34 1 135 1 [58 1 163] 
is provided in Appendix[C]in the case of Example ll.2[ To describe one of the principal differences 
between the two weak maximum principles, we recall the definition of the Fichera function [631 
Equations (1.1.2) &: (1.1.3)] (taking into account our sign convention in ()1.3p for the coefficients 
(a,6, c) of the non-divergence-form operator A)^ 

where n denotes inward-pointing unit normal vector field along dG (now assumed, for example, to 
be C^). By using the Fichera weak maximum principle to decide whether to impose a Dirichlet 
condition for u along S to achieve uniqueness of solutions to (II. ip or ()1.2p . one finds that a 
Dirichlet boundary condition is required for ti on S when 69 < on S, in addition to the usual 
Dirichlet condition on dG\ S, whereas only a Dirichlet boundary condition for u on dG \ S is 
required when 69 ^ on S. 
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However, by instead imposing a Dirichlet boundary condition for u on d^\Ti and a second-order 
boundary condition for u G C'^{ff)nCl-^^{iffUT,)nCioc{^) on S, namely tr{aD'^u) E Ciod^UT.) and 
tr{aD^u) = on S, we achieve uniqueness independent of the sign of the Fichera function, bo, on 
S. Since the second-order boundary condition is automaticahy obeyed by, for example, functions 
in the weighted Holder spaces introduced by Daskalopoulos, Hamilton, and Koch [19\ I47j. we are 
naturally led to a more convenient and powerful framework for establishing existence, uniqueness, 
and regularity of solutions. 

Similar remarks apply in the case of solutions to variational equations and inequalities and the 
framework of weighted Sobolev spaces introduced by Daskalopoulos and the author in [18] and 
by Koch in (47] , 

Stochastic representation (or Feynman-Kac) formulae provide a different approach to estab- 
lishing uniqueness of solutions to (jl.ip or (ll.2p and, in the case of A as in Example 11.21 are 
obtained by Pop and the author in [32j. 

Uniqueness results for solutions to the parabolic porous medium equation (and its linearization) 
were established by Daskalopoulos, Hamilton, Rhee, and Koch in |19( [20| |47] but, unlike the 
linearization of the porous medium, for the coefficients (a, b) in (jl.3p we permit a{x) and x ■ b{x) 
to have quadratic growth with x as x — )• oo and, even when the coefficients 6* are constant, we 
do not require that 6* = when i = 1, . . . ,d — 1. A weak maximum principle for the parabolic 
(model) Kimura diffusion operator is given by Epstein and Mazzeo in \27\ Proposition 4.1.1], 
who also employ a form of second-order boundary condition, together with a Hopf lemma and 
a strong maximum principle in \27\ Lemma 4.2.4 &: 4.2.5]. Related uniqueness results and weak 
maximum principles for classical (sub-)solutions to second-order, linear, degenerate elliptic and 
parabolic operators are proved by Pozio, Punzo, and Tesei in \59 \ \61 \ [62]. but they do not make 
use of second-order boundary conditions. 

A weak maximum principle for variational (sub-)solutions to second-order, nonlinear, degen- 
erate elliptic operators is established by Bonafede [10], with a complex collection of hypotheses, 
while a weak maximum principle for variational (sub-) solutions to (jl.ip is proved by Monticelli 
and Payne |54j when it is assumed that the coefficient matrix, a{x), of A in (jl.3p possesses a 
uniformly elliptic direction and when is bounded (conditions, among others, which we do not 
impose in this article). Borsuk obtains a weak maximum principle, a Hopf lemma, and a strong 
maximum principle |1H Theorem 6.2.1, Lemma 6.2.1 & Theorem 6.2.2] for (sub-)solutions to 
variational equations defined by certain degenerate elliptic quasilinear operators on non-smooth 
domains using weighted Sobolev spaces. Our maximum principle for variational (sub-)solutions 
to (jl.ip (Theorems 18.111 and l8.15"]) generalizes that of Trudinger (TOj Theorem 1]. 

1.5. Notation and conventions. For x,y G M, we denote x A y := uim{x,y} and x \/ y := 
max{x, y}, while x~^ := xV and x~ := —{x AO). We let B{p, R) C M*^ denote the open ball with 
radius R > and center p G M"^. For an integer k > 0, we let C^{i?) denote the vector space of 
functions whose derivatives up to order k are continuous on 0' and let C^{ff) denote the Banach 
space of functions whose derivatives up to order k are uniformly continuous and hounded on ^ [3l 
§1.25 & §1.26]. If T ^ di? is a relatively open set, we let C^^d^ ^ denote the vector space of 
functions, n, such that, for any precompact open subset C/ d U T, we have u G C^{U). When 
we label a condition an Assumption, then it is considered to be universal and in effect throughout 
this article and so not referenced explicitly in theorem and similar statements; when we label a 
condition a Hypothesis, then it is only considered to be in effect when explicitly referenced. 

1.6. Acknowledgments. I am grateful to Panagiota Daskalopoulos, Peter Laurence, and Camelia 
Pop for many engaging discussions on degenerate partial differential equations. 
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Part 1. Weak and strong maximum principles for operators on smooth functions 
and applications to boundary value and obstacle problems 

In this part of our article 33 SI 35]), we develop weak and strong maximum principles 
for operators on smooth functions and their applications to boundary value and obstacle problems. 

2. Applications of the weak maximum principle property to boundary value 

problems 

We shall encounter many different situations (for example, bounded or unbounded domains 
& <ZW^, bounded or unbounded functions u with prescribed growth, and so on) where a basic 
maximum principle holds for linear, second-order, partial differential operators A in (|1.3|) with 
non-negative characteristic form and acting on a linear subspace of functions in C^(^). In 
order to unify our treatment of applications, we find it useful to isolate a key "weak maximum 
principle property" (Definition I2.8P and then derive the consequences which necessarily follow 
in an essentially formal manner. In this section we consider applications to boundary value 
problems. After providing some technical conditions in §2.11 we proceed to the main applications 
in §2.21 namely a comparison principle for subsolutions and supersolutions and uniqueness 
for solutions to the Dirichlet boundary problem (Proposition I2.16P and a priori estimates for 
subsolutions, supersolutions, and solutions (Proposition 12.19]) . In §2.31 we show that when 
operator has the weak maximum principle property for subsolutions which are bounded above, 
the property may also hold subsolutions which instead obey a growth condition (Theorem 12. 20p . 

2.1. Preliminaries. We introduce the 

Definition 2.1 (Non-negative definite characteristic form and boundary degeneracy locus). Sup- 
pose a = (a*{} is an R"^^ '^-valued function on d' which defines a non-negative definite characteristic 
form, that iqj 

{ar],r])>0 on ^, Vr? G M"'. (2.1) 

We call 

j::=mt{{xedff:a{x)=0})gdff (2.2) 
the boundary degeneracy locus defined by the coefficient matrix a = (a*-'). □ 

In this article, we normally assume that S in (|2.2p is non-empty. Of course, when the matrix 
a{x) in Definition 12.11 is symmetric for all x G dG , then the definition of S in ()2.2p is equivalent 
to the definition 

■= {xGdff : {a{x)r],r]) = 0,Vr? € M'^}, 

used by Fichera [35j . Olemik, and Radkevic [58], [631 P- 308] (where the whole domain boundary 
is denoted by S). As in [63] §1-1] , we have the 

Definition 2.2 (Linear, second-order, partial differential operator with non-negative definite 
characteristic form). We call A in ()1.3p a linear, second-order, partial differential operator with 
non-negative definite characteristic form on G if the coefficient matrix, a, obeys ()2.ip . 



We follow the definition of Olemik, and Radkevic |63l p. 308] rather than Tricomi |63l p. 298], which requires 
in addition that (a?7, ??) > on i^t" for all ?7 G K'* \ {0}; however, in the applications we have in mind, the latter 
condition is always satisfied and the degeneracy is confined to a subset of dG . 



PARTIAL DIFFERENTIAL OPERATORS WITH NON-NEGATIVE CHARACTERISTIC FORM 



11 



Definition 2.3 (Second-order boundary condition). Let a: G ^ M'^^'^ obejj^ 

a G Cioc(^US;M'^^'^). (2.3) 

We say that u G C'^{ff) R Cioc(^) obeys a second- order boundary condition along S defined by 
a G aoc(^;M'^'"^) via ([22]) if 

L>n G Cioc(^US;M^), (2.4) 
tr(aL>2u) G Cioc(^ U S), (2.5) 
iT{aD'^u) = on S. (2.6) 

Remark 2.4 (Second-order boundary conditions and Holder spaces). For certain second-order 
differential operators with non-negative characteristic form, such as the linearization of the porous 
medium equation, the second-order boundary condition (|2.6|) is a property of functions in suitable 
Holder spaces, C|^"(^ U S), defined in [19] for functions on a domain ^ <zW^ equipped with a 
Riemannian metric g^^ = a*-' , where a is as in Definition 12.11 and we replace the usual Euclidean 
distance function on in the definition of standard Holder spaces by the Riemannian distance 
function, s, defined by the metric a^K See [IHl Proposition 1.12.1], [301 Lemma 3.1] and [29] for 
further discussion. □ 

Remark 2.5 (Ventcel boundary conditions and Feller semigroups). The second-order boundary 
condition in Definition 12.31 may be viewed as a special case of a generalized Ventcel boundary 
condition |681 §7.1]. Around 1950, Feller completely characterized the analytic structure of one- 
dimensional diffusion processes, giving an intrinsic representation of the infinitesimal generator 
^ = —A of a one-dimensional diffusion process and determined all possible boundary conditions 
which describe the domain, ^{s^). (A parallel analysis from the point of view of Sturm-Liouville 
operators may be found in [72].) In 1959, Ventcel studied the problem of determining all possible 
boundary conditions for multi-dimensional diffusion processes and a generalization of Ventcel's 
results to Feller processes, by Taira and many others, is described in [681 Chapter 7]. However, 
while the Ventcel boundary condition may be exogenous (that is, defined by an auxiliary, degen- 
erate second-order operator, L, which is distinct from A), the second-order boundary condition 
in Definition 12.31 are endogenous (that is, naturally induced by j4). □ 

Remark 2.6 (Relationship with definition of degenerate elliptic operators in the theory of vis- 
cosity solutions). Consider a second-order partial differential equation of the form 

F{x, u, Du, D'^u) = 0, Vx G ^, 

defined by a map 

F : xRxR'^ X S'^ ^R, 

where C M"'^'^ denotes the subset of real, symmetric matrices. The map, F, is called degenerate 
elliptic (in the sense of viscosity solutions) [15\ p. 1] if it is increasing with respect to its matrix 
argument, 

F{x, r, 7], X) < F{x, r, r],Y), \/X <Y and (x, r, r/) G ^ x M x M'^, 
where X < y in the sense of S"^ if 



"^It would be more than sufficient for tlie application of Definition 12.31 in this article to require only that a is 
continuous on some open neighborhood of E in ff. 
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When F is degenerate elliptic in the preceding sense, the map is called proper if 

F{x, r, r], X) < F{x, s, rj, Y), yX <Y,r<s and [x, ??) E ^ x R"'. 

Of course, the linear operators A in (II. 3|) and equations considered in this article will always 
define degenerate, proper maps F{x, u, Du, D^u) := Au — f since the coefficient matrix a of D^u 
will obey (12. ip and the coefficient of u will obey c > on however, the converse is not true 
— when F is degenerate elliptic in the sense of [151 p. 1], it does not imply that A is degenerate 
elliptic in the sense of this article (meaning a lack of strict inequality in (II. 3p on ^). □ 

To simplify the statements of our definitions and results, we shall make the universal 

Assumption 2.7 (Partial differential operator with non-negative characteristic form). Given 
a domain, 0' ^ M"^, the coefficients of the linear, second-order partial differential operator A : 
C^{ff) C(^) as in ([O]) obey ([21]), ([231), and the open subset, S g d^, is given by dO]) . 

We can now state the key 

Definition 2.8 (Weak maximum principle property for subsolutions). We say that an oper- 
ator A in (II. 3p obeys the weak maximum principle property on ^ U S if whenever u G C^(i^) H 
Cioc(^) obeys ([23]), dM]) anci 

Au<0 on ^, (2.7) 
n<0 on(9^\S, (2.8) 
and, if & is unbounded, sup^u < oo, then 

U<0 OVl d. 

Examples of operators, A^ with the weak maximum principle property on ^ U S are provided 
by Theorems EH [531 EiM and (when S = 0) gOl Theorem 3.1 & Corollary 3.2]. 

Remark 2.9 (Auxiliary boundedness condition in the definition of the weak maximum principle 
property). Once A is known to have the weak maximum principle property on i^US, the condition 
that u be bounded above when G is unbounded (which could be replaced by a growth condition 
via Theorem l2.20p plays no role in most applications of the weak maximum principle property. □ 

Remark 2.10 (Auxiliary conditions required for linear, second-order elliptic differential operators 
to have the weak maximum principle property on S). In order to prove that A has the weak 
maximum principle property on U S, we shall need to define open subset, S ^ 36 ^ via ()2.2p 
and require that A to obey auxiliary conditions such as ()2.ip . ()2.14p . (j2.13p . ()5.2p or (j5.3p . ()5.4p . 
()5.7p in order to prove that A has the weak maximum principle property on €?US. However, once 
A is known to have the weak maximum principle property on i^U S, neither the definition of the 
open subset, S ^ dG ^ via ()2.2p nor these auxiliary conditions play any role in most applications 
of the weak maximum principle property. □ 

Remark 2.11 (Role of the second-order boundary condition). We chose to make the auxiliary 
conditions ()2.4p . (j2.5p . ()2.6p for u along S part of the definition of the function space for u and, 
as in |19t I30j . they are properties of functions in C|''~"(^ US). In the application of the weak 
maximum principle property, the second-order boundary condition plays no explicit role but we 
implicitly need to restrict functions u G C^(^) n Cioc(i^) to some linear subspace where the 
replacement of dG by dG \ S makes sense and the subclass of functions obeying a second-order 

^The condition (|2.6p could be replaced by tr(a_D^u) < on E, for a slight increase in generality, as is clear from 
Step [T] of the proof of Theorem 15.11 
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boundary condition provides the arena in which we shah give examples of operators, A, having 
the weak maximum principle property. □ 

Remark 2.12 (Non-negative zeroth-order coefficient). The condition ()2.13p on the coefficient of 
u, namely, 

c > on 

is not explicitly required in Definition 12 . 81 since A is already assumed to possess the weak maximum 
principle property on S, although simple counterexamples to the weak maximum principle 
exist when this condition is relaxed in general \4:0\ p. 33]; the condition c > on ^ is a necessary 
condition for the strong maximum principle [60^ Exercise 2.1]. □ 

Remark 2.13 (Equivalent first-order boundary conditions). If n G C'^{ff)riCl^^{^UT,)nCioc{^) 
obeys (12. Sp and (12. 7p . with A as in in (|1.3p . then the second-order boundary condition (12. 6|) is 
equivalent to 

- {b, Du) + cn < on S. (2.9) 
Indeed, when we have equality in (j2.7p . and thus equality in (j2.9p . the condition (|2.9p is analo- 
gous to the boundary condition proposed by Heston [ISj Equation (9)] for the parabolic termi- 
nal/boundary problem corresponding to the boundary value problem (jl.ip : one obtains 

-{b,Du)+cu = f onS, (2.10) 

for the boundary value problem (jl.ip when / is non-zero and u G C^(^) H Cjq^(^U S) n Cioc(^) 
obeys (12. 5p . Indeed, the parabolic version of the condition (j2.10p (normally when / = 0) is often 
used in the numerical solution of parabolic boundary value or obstacle problems in mathematical 
finance [Ml Equation (22.19)], [731 Equation (15)]. □ 

2.2. Applications of the weak maximum principle property to boundary value prob- 
lems. We begin with the 

Definition 2.14 (Classical solution, subsolution, and supersolution to a boundary value problem 
with Dirichlet data). Given functions / G C(^) and g G Cioc{d^\^), we call u G C^(^)nCioc(^) 
a subsolution to a boundary value problem for an operator A in (II. 3p with Dirichlet boundary 
condition along dff \ S, if u obeys ()2.4p . ()2.5|) . (12. 6p . sup^u < oo, and 

Au< f on (2.11) 

u<g on^\S. (2.12) 

We call u a supersolution if —u is a subsolution and we call u a solution if it is both a subsolution 
and supersolution. □ 

When S = 0, then Definition 12 . 141 reverts to the standard definition of a solution to a boundary 
value problem [50]. The first application, of course, of the weak maximum principle property is 
to settle the question of uniqueness for solutions to the Dirichlet boundary problem. 

Remark 2.15 (Auxiliary boundedness condition in the definition of subsolution, supersolution, 
and solution). The condition in Definition 12.141 that a subsolution, u, be bounded above when 
^ is unbounded can be replaced by a growth condition via Theorem 12.201 and similarly for a 
supersolution or solution. □ 

Proposition 2.16 (Comparison principle for subsolutions and supersolutions and uniqueness for 
solutions to the Dirichlet boundary problem). Let A in (jl.3p obey the weak maximum principle 
property on ffUT,. Let f G C(^) and g G C\oc{dff \ S). Suppose ui is a subsolution and U2 a 
supersolution to the Dirichlet boundary problem in the sense of Definition \2.14\ Then U2 > ui 
on G and, if ui,U2 are solutions, then U2 = ui on G . 
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Proof. Since ui is a subsolution and U2 a supersolution, then ui — «2 is a subsolution with 
«i — it2 = on dff \ S and thus 

ui — U2 <0 on S', 

because A has weak maximum principle property on ^US. When ui,U2 are both solutions, then 
we also obtain U2 — ui < on 0' and U2 = ui on 0". □ 

We shall usually need to appeal to one of the following conditions. 

Hypothesis 2.17 (Non-negative zeroth-order coefficient). The coefficient c in ()1.3p obeys 

c > on ff. (2.13) 

Hypothesis 2.18 (Uniform positive lower bound on zeroth-order coefficient). There is a constant 
Co > such that the coefficient c in (II. Sp obeys 



c > Co on ff. (2.14) 

We can now proceed to give the expected a priori estimates. 

Proposition 2.19 (A priori estimates for subsolutions, supersolutions, and solutions). Let 
A in (jl.3p obey the weak maximum principle property on ^UTi and assume c obeys (|2.13p . Let 
f G C(i^) and g £ Ciodd^X S). Suppose u is a subsolution to the boundary value problem in the 
sense of Definition \2.14\ 

(1) // / < on and u is a subsolution for f and g, then 

n < V sup g on G . 

(2) If f has arbitrary sign and u is a subsolution for f and g but, in addition, there is a 
constant co > such that c obeys (I2.14p . then 

ti < V — sup / V sup g on d . 

Co e diff\T. 

(3) If f > on 0" and u is a supersolution for f and g, then 

u > A inf g on G . 

(4) If f has arbitrary sign, u is a supersolution for f and g (Definition \2. and c obeys 
(HIH), then 

u > A — inf f A inf g on G. 

Co ff d&\T, 

(5) If f = on G and u is a solution for f and g (Definition \2. 14^ , then 

\\u\\c(ff) ^ \\9\\c(de\E)- 

(6) If f has arbitrary sign, u is a solution for f and g, and c obeys (j2.14p . then 

\\'^\\c{ff) ^ —\\f\\c{d-) V Il5'llc{a^)' 

CO 

The terms supg^\^2 9? '^^'^ iiifai^\s9; ^^'^ ll^llc( (9^\s ) preceding items are omitted when 
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The a priori estimate in Item ([6]) may be compared with \63\ Theorem 1.1.2] (in the case of 
functions) and [63^ Theorem 1.5.1 & 1.5.5] and [69^ Lemma 2.8] (in the case of functions). 
However, because the coefficient matrix a of A in (jl.3p is zero along S ^ dff in apphcations 
considered in this article, there is no analogue of [40 ^ Theorem 3.7 & Corollary 3.8]. 

Proof of Proposition \2.1fA For Items ([T]) and ([2]), we describe the proof when S ^ dff] the proof 
for the case S = di? is the same except that the suprema on the right-hand side are replaced by 
zero. When / has arbitrary sign, choose 

M := V — sup / V sup g, 

Co ff d£i\T, 

while if / < 0, choose 

M := V sup g. 

We may assume without loss of generality that M < oo. Since M > g on dff \ Ti, we have 

AM = cM>coM>f on ^ (by (EH])), 

when / has arbitrary sign and, when / < 0, we have 

AM = cM > > / on ^ (by (l2l3D ). 

Thus, M is a supersolution and so u < M on ^ by Proposition 12.161 which gives Items ([T]) and 
([2]). Items (l3|) and follow from Items ([T|) and ^ by noting that — n is a subsolution for — / 
and —g if u is a supersolution for / and g. Item ([5]) follows by combining Items ([1]) and ([3]), while 
Item dni) follows by combining Items ([2]) and dH). □ 

2.3. Applications of the weak maximum principle property to unbounded functions. 

Our Definition l2.8l of the weak maximum principle property requires that the subsolution, u, on an 
unbounded domain is bounded above. However, if an operator has the weak maximum principle 
property for subsolutions which are bounded above, we obtain an extension for subsolutions which 
instead obey a growth condition. 

Theorem 2.20 (Weak maximum principle for unbounded functions on unbounded domains). 
Let A be an operator as in (II. 3p . let ^ M'^ 6e a possibly unbounded domain with T, ^ as in 
([22]), and let (p G C'^{ff) obey < tp < 1 on ^ . Let 

Bv := -[A,ip]{ip~^v), yv£C'^{ff), (2.15) 

and suppose that the differential operator, 

A:={A + B)v, yveC^iff), (2.16) 

has the weak maximum principle property on ffUT, (Definition \2.8\) for functions u on which 
are bounded above. Then A has the weak maximum principle property on ffUT, for functions u 
on G which obey the growth condition, 

u<C{l + ip-'^) onG. (2.17) 

Proof. Clearly, we have 

A{^pu) = Aipu + Bifu 

= ipAu + [A, ip]u + B(pu 

< on ^ (by dZZ]) and dUSD). 
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Since ipu < C{ip + 1) < 2C on i^, then (pu < on 0' since A has the weak maximum principle 
property on U S for functions u on ^ which are bounded above. Thus, u < on 0'. □ 

3. Applications of the weak maximum principle property to obstacle problems 

We now turn to the apphcation of the weak maximum principle property to obstacle problems. 
The application is complicated by the fact that the optimal interior regularity of solutions to ob- 
stacle problems is C^'^(^) rather than C'^(ff). For this reason, we first consider in ^3.1l a simple, 
if artificial, special case where the solutions are assumed to be C^, and prove a comparison prin- 
ciple for supersolutions and uniqueness for solutions to the obstacle problem (Proposition 
13. Sp and derive a priori estimates for supersolutions and solutions to obstacle problems 
(Proposition 13. 5( ). 

In ^3.2| we develop a theory of continuous supersolutions to the obstacle problem, modeled 
on the corresponding theory for the Dirichlet boundary value problem in |40[ §2.8 &: §6.3] for 
linear, second-order, uniformly elliptic partial differential equations on bounded domains. Our 
development has some independent interest since it is a theory of classical rather than viscosity 
solutions, unlike that of |16] . We then consider the full problem of developing applications of 
the weak maximum principle property to C^'^(^) solutions to obstacle problems. In particular, 
we prove a comparison principle for continuous supersolutions and uniqueness for solutions to 
the obstacle problem (Theorem I3.15P and a weak maximum principle and a priori estimates for 
solutions and continuous supersolutions to the obstacle problem (Proposition 13. 18]) . 

3.1. Applications of the weak maximum principle property to obstacle problems with 
supersolutions. We have the 

Definition 3.1 (Solution and supersolution to an obstacle problem). Given functions / E 
C(^), and g £ aoc(5^\S), and ^ G Ci'^(^) n Cioc(^) obeying V < 5 on dff\E and (ffs]) . 
(12. 6p . we call u £ C^'^(^) R Cioc(^) a solution to an obstacle problem for an operator A in (II. 3p 
with Dirichlet boundary condition along dff \ S, if n obeys (12. 4p . (j2.5p . (j2.6p . sup^ |n| < oo, and 

nim{Au — f,u — Tp} = a.e. on ff, (3.1) 
u = g on^\S. (3.2) 

We call n G C2(^) n Ciod^) a supersolution to the obstacle problem if u obeys ([Ti|) . (j23|) . 
()2.6p . supffU > — oo, and ()3.ip and ()3.2p hold with = replaced by >. 

The comments in Remark 12. 151 applv equally well to Definition 13.11 

Remark 3.2 (Regularity properties of the obstacle function). We require the obstacle, ip, in 
Definition 13. II to share the same regularity properties as those of the solution, u, since (unless the 
obstacle problem is trivial, in the sense that u> ip solves Au = f on ff), the subset 0'r]{u = ip} 
will be non-empty. However, there are many examples in applications where ip S Cj^'^ {^) but 
ip = ipi f\'il)2 and 1^1,1^2 G C'j^^(i^) are convex (such as the perpetual American put [Ml §8.3]) and 
the obstacle problem in Definition 13. II can still be well-posed in such cases. 

Proposition 13.31 below for solutions to the obstacle problem is an analogue of Theorem 17. 2( 
which applies to i/^(^, tt)) solutions to a variational inequality, and Theorem 17.71 which applies 
to H'^{ff, tv), that is, strong solutions to the obstacle problem. We may also compare Propositions 
ED and ESI with [Ml Theorems 4.5.1, 4.6.1, 4.6.6, & 4.7.4, and Corollary 4.5.2] for the case of 
variational inequalities. 
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Proposition 3.3 (Comparison principle for supersolutions and uniqueness for solutions 
to the obstacle problem). Let A in (jl.3p obey the weak maximum principle property on S. 
Let f G C{ff), and g £ Cxoc{d& \ S), and ^ G C^{ff) n QocC^) with i!<gondff\^. Suppose 
ui is a solution and U2 is a supersolution to the obstacle problem in the sense of Definition 
\3.1\. Then U2 > ui on d' and if U2 is also a solution, then U2 = ui on G . 

Proof. Suppose ui is a solution and U2 a supersolution to the obstacle problem and that := 
n {ui > U2} is non-empty. Observe that 5^ = {^n{ui = U2}) U {dGn{ui > U2}) and 
because ui = g < U2 on d& \ S, we must have d& H {ui > U2} = S H {ui > M2}. Therefore, 
ui — U2 < on 9^ \ S. We have ui G C^(^) by elliptic regularity [Ml Theorem 6.13] and 
U2 G C^('^) because U2 is a supersolution, so ui — U2 £ C'^{'W) and ui — U2 obeys (I2.4p . (j2.5p . 
(I2.6P with in place of ff. Moreover, A{ui — U2) < on and so ui — «2 < on by the 
weak maximum principle property for A (Definition 12 .Sp . contradicting our assertion that ^ is 
non-empty. Hence, ui < U2 on 0" and thus ui < U2 on G . □ 

Remark 3.4 (Relaxation of the C^(^) regularity condition in Propositions 13.3] and [33]) . We 
include Propositions 13.3 1 and 13.51 because their proofs short, transparent, and insightful, but a far 
more useful versions are provided by Propositions 13.15) and [3T8l respectively, where the C^(^) 
restriction is removed. 

We then have the 

Proposition 3.5 (Weak maximum principle and a priori estimates for supersolutions and 
solutions to obstacle problems). Let A in (|1.3p obey the weak maximum principle property 
on G U and assume c obeys (j2.13p . Let f G C{G), and g G CioddG \ S), and G Cioc(^) 
with tp < g on 36 \ S. Suppose u is a C'^ supersolution to the obstacle problem in the sense of 
Definition \3.1[ 

(1) lff>0 on 6, then 

u> A inf g on & . 

(2) If f has arbitrary sign but there is a constant cq > such that c obeys (I2.14p . then 

u > A — inf / A inf o on ff. 

Co d6\^ 

(3) If f < on and u is a C"^ solution for f and g and ip, then 

n < V sup g V sup ip on d . 

(4) If f has arbitrary sign, u is a solution for f and g and ijj, and c obeys ()2.14p . then 

u < V — sup / V sup g V sup -0 on G . 

(5) // ui and U2 are solutions, respectively, for fi > f2 and ipi > 'ip2 on G , and gi > 52 
on 80 \ S, then 

ui > U2 on ff. 

(6) If ui and U2 are C"^ solutions, respectively, for fi,ipi and f 2,1^2 on G, and gi and g2 on 
80' \T., and c obeys ([ZTip . then 

hi - Mciff) ^ - /2|lc(^^) V llffi -ff2|lc(a^) V IIV'i -V'2|lc(<f)> 
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while if fi = / = /2 o,nd c obeys (|2.13p . then 

\\ui - U2\\c(ff) < \\gi - 92\\c(aws^) V IIV'i - V'2|lc(^)- 

Proof. Consider Items ([T]) and Since u is a supersolution to the obstacle problem in Definition 
I3.H then it is also a supersolution to the boundary value problem in Definition 12.141 (where "0 
plays no role) and so Items (JT]) and ([2]) here just restate Items ([3]) and (jl} in Proposition 12.191 
Consider Items ^ and ^ here. When / < on and ^ obeys (12131) . let 

M := V sup g V sup ijj, 

while if / has arbitrary sign and A obeys (I2.14p . let 

M := V — sup / V sup g V sup 

Co d0\T. e 

We may assume without loss of generality that M < oo. Then M > -0 on ^ and M > (7 on 
S, while 

AM = cM >0> f on ff, 

when / < on ^ and 

AM = cM > cqM > f on ^, 

when / has arbitrary sign. Hence, M is a supersolution and so Proposition 13 . 31 implies that u < M 
on 1^, which establishes Items ^ and For Item ([5]), observe that ui is a supersolution for 
the obstacle problem in Definition 13.11 given by f2,g2,tp2 and thus ui > U2 on ^ by Proposition 
[331 For Item §1, define 

:= —ll/i - /2|lc{<?) V llfi-i - 5211^(9^) V IIV'i - V'2|lc(^?) and u := U2 + m. 

Then 

Au = Au2 + Am>f2 + cm>f2 + sup(/i - /2) > /i on ff, 

Iff 

while 

u > ^2 + sup(V'i — ^2) > V'l oil '^j 

and 

u>g2+ sup (51 - 52) > 51 on c^i^ \ S. 

Therefore, w is a supersolution for fi,gi,ipi and so Proposition 13.31 implies that u > ui on ff, 
and thus 

Ui — U2 < m on G . 

By interchanging the roles of tti, U2 in the preceding argument, the conclusion follows for the case 
c > Co > 0. For the case c > 0, we now define 

"T- := II V'l - V'2|lc(^) V llfi-i -5211(7(9^) and u := ^2 + m, 

so that 

Au = Au2 + Am > / + cm > f on 0', 
and the remainder of the argument is identical. □ 
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3.2. Applications of the weak maximum principle property to obstacle problems with 
continuous supersolutions. Before introducing our definition of a continuous supersolution, 
we shall need the 

Definition 3.6 (Local models for open neighborhoods of points in ^). Suppose that G C M"^ is 
a domain with boundary, that 30 is a submanifold-without-boundary in M"^ of dimension (i, 
and that S C dG is a submanifold in M*^ of dimension d — 1 with boundary 9S being a 
submanifold-without-boundary in M"^ of dimension d—2. We call C/ d S a model neighborhood 
of a point p G U S if there is a C^-diffeomorphism $ : M'^ — ^ M'^ such that ^{p) = and 

(1) If p G 1^, then #([/) = 5(0, r), where 

5(0, r) := {x G M'^ : |x| < r} 

is the open ball of radius r > 0; 

(2) If p G S, then $([/) = 5+(0,r), where 

5+(0,r) := {x G M'' : |x| < r,Xrf > 0} C W^~^ x M+ 
is the open half-ball, and 

$(C7 n S) = 5(0, r) n {xd = 0}; 

(3) If p G as, then ^>(C/) = 5++(0,r), where 

5++(0,r) := {x G M'^ : \x\ < r,Xd-i > 0,Xrf > 0} C R''"^ x R+ x R+ 
is the open quarter-ball, and 

$(t7n(a^\S)) 
$(i7ns) 
$(^n5S) 

We shall also need the 

Hypothesis 3.7 (Existence and uniqueness of solutions to the Dirichlet problem on model neigh- 
borhoods). Suppose that [/ C R'^ is a model neighborhood in the sense of Definition 13. 9| that / 
and the coefficients of A in ()1.3p belong to C^{U) for some a G (0, 1), that S = {x G dU : a{x) = 
0}, that g G C{dU \ S), and that A has the weak maximum principle property on U S in the 
sense of Definition 12.81 

Then there is a unique solution u G C'^~^°'{U) H C{U) to the Dirichlet problem, 

Au = f on U, 
u = g on dU \ S, 

which obeys ([MD, (l23|), and ([M]). 

Example 3.8 (Examples where the conclusions of Hypothesis 13.71 hold). When U is a ball and 
T, = and c > on ?7, then |40l Theorem 6.13] yields the conclusion of Hypothesis 13.71 for 
the resulting uniformly elliptic operator A on U. Although not a model neighborhood, since 
U in Hypothesis 13.71 is assumed bounded, the articles [HI [20l [30] provide examples of parabolic 
differential operators with non-negative characteristic form where the conclusion of Hypothesis 
[321 holds for U = R'^"'^ x R+, and S = R^^^^ x {0}, and b'^ > 0. 



= 5(0,r)n{xrf_i = 0}n{xrf>0}, 
= 5(0,r)n{xrf„i >0}n{xd = 0}, 
= B(0,r)n{xd-i=Xd = 0}. 
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By analogy with the definitions of continuous superharmonic functions [401 §2.8], [441 Definition 
6.1] or continuous supersolutions to linear, second-order elliptic partial differential equations [iQl 
pp. 102-103], we have thd 

Definition 3.9 (Supersolution to an obstacle problem). Let ^ and S ^ 9^ be as in 
Definition 13.61 and let A in (jl.3p obey Hypothesis 13.71 Given functions / € C°(^), and g € 
C\oc{d& \ S), and £ C'ioc(^) obeying ip < g on dff \ S, we call v € C(^) a supersolution to 
the obstacle problem (j3.ip . (j3.2p if u > -0 on ^ and inf^f > — oo and if for every open subset 
[/ (E U S) n {u > V'} such that U is model neighborhood in the sense of Definition 13.61 and for 
every v G C"^ {U) r\ C {tj) obeying 

Av = f on U, 
V > V on dU \ E, 

and ([23]), (US]), (^M), we then have 

V > V on U. 

We call V G Cioc(^) a supersolution relative to g on dG \ S if v is a supersolution and v g on 

Clearly, every solution or C'^{tff) n Cioc(^) supersolution to the obstacle problem in the sense 
of Definition 13. II is also a Cioc(^) supersolution relative to g. Note that v in Definition 13.91 is not 
required to obey v > "0 on ?7 and that, unlike in Definition 13. 11 we only require that the obstacle, 
■0, be continuous on ^. 

Remark 3.10 (Continuous supersolution and subsolution to a boundary value problem). When 
the obstacle, is omitted (in the obvious way) from Definition 13.91 then we call v a supersolution 
(relative to g on dG \ Ej to the boundary value problem (jl.ip . Similarly, we call v a subsolution 
(relative to g on di? \ to the boundary value problem (11. ip if —v is a supersolution (relative 
to g on S). 

Remark 3.11 (Perron's method and solution to an obstacle or boundary value problem). It 
would be natural to seek a solution to the obstacle problem in Definition 13.11 via Perron's method 
(see [lOl Theorems 2.12 & 6.11]) and attempt to show that 

u{x) = inf v{x), a; G ^, (3.3) 

is a solution to the obstacle problem, where '^f^g^ip is the set of all supersolutions to the obstacle 
problem in Definition 13.91 defined by /, 5, and ip. We shall not consider the question of existence 
itself in this article, beyond providing an analogu^ of [44^ Lemma 6.3] in Lemma 13.12^ but in 
Theorem l3.151 we show that if « is a solution then it can be expressed in the form (13. 3p . Naturally, 
this method would also apply to the boundary value problem (jl.lj) . 

Lemma 3.12 (Construction of smaller supersolutions to obstacle problems). Let G and 

Ti ^ diT" he as in Definition \3.(^ and let A obey Hypothesis \3.'T\ Let f £ C°(^), and g G 
Cioc{dff \ S), and tp G Cioc(^) obey ip < g on di^ \ S. Suppose that v G C(^) is a supersolution 
to the obstacle problem in the sense of Definition \3.9l Let U (£ {ff L) 'S) Ci {v > ip} be a model 



''While we could define supersolutions to the obstacle problem in the viscosity sense — see [161 Example 1.7 & 
Definition 2.2] and [51| — this simpler and more tractable definition will suffice for our applications. 

^Lemma 6.3 in [33] asserts that if t; £ C{ff) is subharmonic, B <£ ff, and w is defined by Aw = on B and 
w = V on ^ \ B, then w is subharmonic and v < w on ff. 
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neighborhood in the sense of Definition \3.6l Define v G C{^) by setting v = v on 0' \ U and 
choosing v G C^+°(C/) H C{tJ) to be the unique solution to 

{Av = f on U, 
V = V on dU \ S. 

// at least one of the following conditions hold, 

(1) c> on U and f > on U and ip >0 on U ; 

(2) c> cq > on U for some positive constant cq; 

then V is a supersolution for small enough r > depending on p £ S and tp. 

Remark 3.13 (Application to Dirichlet boundary value problems). Lemma 13.121 also yields 
supersolutions to the Dirichlet boundary value problem (ll.ip . 

Proof. Consider ([TJ. Since / > on [/, then Av > and the weak maximum principle property 
implies that 

> A inf v. 

dU 

Since v is continuous on 0" and v{p) > il^{p), we may choose r in Definition 13.61 small enough that 

inf v > sup ip. 

dU u 

Since ip > on U and thus v > on U, we obtain 

inf V > inf v > sup ih, 

u au u 

and hence 

V > on U. 

But v = v>'iIj on ff\U and therefore v is a supersolution. 

Consider ([2]). In this case, / has arbitrary sign but c > cq > and so we may choose a 
positive constant m such that / + cm > f + cqth > on C/ and ip + m > on [/. Then 
A[v + m) = Av + cm = / + cm >Q onU and the conclusion for Item ([T|) yields v + m > + m 
on U and thus we again obtain v > ^p on U, as desired. □ 



We introduce an analogue of Definition [27 

Definition 3.14 (Strong maximum principle property for functions). Let C M'' be a domain 
and let S Q dG be as in ()2.2I) . We say that an operator A in (jl.Sh obeys the strong maximum 
principle property on ^ U S if whenever u G C^(^) n Cioc(^) obeys (|2.4p . (12.5P and (I2.6j) and 
Au < on 0', then one the following holds. 

(1) If c = and u achieves a maximum in U S, then u is constant on ff. 

(2) If c > and u achieves a non-negative maximum in S, then u is constant on S". 

Theorem 14.61 and (when S = 0) [401 Theorem 3.5] give examples of operators, A, with the 
strong maximum principle property. 

Theorem 3.15 (Comparison principle for supersolutions and uniqueness for solutions to the 
obstacle problem). Let ^ ^ M'^ and S ^ dff be as in Definition \3.6\ and let A in ()1.3p obey 
Hypothesis \ 3.1\ and have the strong maximum principle property on S. Let f G C"(^), and 
g G C\ac{di^ \ Ti), and ip G Cioc(^) obey ip < g on dG \ S. Suppose that u is a solution to the 
obstacle problem in the sense of Definition \3.1\ and v a supersolution to the obstacle problem in 
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the sense of Definition \3.A In addition, assume that each connected component C of ^n{u > ip} 
has the property that 

9Cn(5^\Sn{5>^}U{n = V'}) (3.4) 
Then u < v on if , and if v is also a solution, then u = v on G . 

We immediately obtain the following analogue of Proposition 12.16] where the hypothesis (13. 4p 
is replaced by the simpler condition that S ^ dG and noting that the proof of Proposition 12.161 
easily adapts to the case where u is a subsolution in the sense of Remark 13.101 rather than a 
solution in the sense of Definition 12.141 

Corollary 3.16 (Comparison principle for continuous supersolutions and uniqueness for solutions 
to the Dirichlet boundary value problem). Assume the hypotheses of Theorem \3.15\ on G , S, A, 
f , and g. Suppose that u is a subsolution and v a supersolution to the boundary value problem 
(jl.ip in the sense of Remark \3.10[ . In addition, assume that 

dG\T.^2S. (3.5) 

Then u <v on G, and ifu,v are both solutions, then u = v on G. 

Proof of Theorem \3.15\ We may assume without loss of generality that 0'r\{u > tj:} \s non-empty, 
as otherwise we are done. Observe that u = g < v on dG \ S and so n — u < on dG \ S, noting 
that, by hypothesis, dG \ S is non-empty. Let 

M := sup(ti — v), 

noting that M < oo by hypothesis (Definitions 13.11 and 13. 9p , and consider 

D := {G [JT.) r\ {u > ^} r\ {u - V = M}. 

Clearly, D is a relatively closed subset of {GVJ S) n {u > ^} in since u, v, and are continuous 
on G. We first assume that (i^ U S) n {u > V} is connected and aim to show that D is also a 
relatively open subset of {G UT.) Ci {u > 

We now proceed by adapting the proou of [44 ^ Lemma 6.2]. Let p ^ GuT, and U d (^U S) n 
{u > "0} be a model neighborhood in the sense of Definition 13.61 By Hypothesis 13.71 there is a 
unique solution v G C2+"([/) n C{U) to 

Av = f on U, 
V = V on dU \ S. 

Then v < v on U (recall that v is a supersolution) and hence 

u — V > u — V on U. 

Elliptic regularity, via Hypothesis [3771 ensures that u G C^+°(C/)nC(C7) since U d {GuT,)n{u > 
ip} and thus Au = f on U. We have A(u — v) = Au — Av = on U and u — v = u — v<Mon 
dU, that is, 

jA{u-v) = onU, 

\u-v<M on dU\T.. 

By the weak maximum principle property for A and the fact that M > 0, we must therefore have 
u — V < supgij\-£{u — v)~^ < M on [/; in particular, 

M > {u- v){p) >{u- v){p) = M. 

^Lemma 6.2 in [S] asserts that if ff C R'* is a bounded domain, u £ C{ff) is subharmonic, v £ C{ff) is 
superharmonic, and u < w on dff, then u <v on ff. 



PARTIAL DIFFERENTIAL OPERATORS WITH NON-NEGATIVE CHARACTERISTIC FORM 23 

Hence, {u — v){p) = M and so u — v has a maximum at p £ C/US. The strong maximum principle 
property for A implies that u — v = M (constant) on U and thus U <Z D. 

Since (^n{ti > ^} is connected by our assumption and D \s a, relatively open and closed subset, 
then either D = or D = n {ti > ^}. If is empty, then u < w on n {u > V'} and u < v 
on n {u = ip}, so n < f on and thus u < v on 0'. 

If D = 0" n {u > ip}, then u — v is constant on D = Ci {u > tp} and thus u — v is constant 
on n {u > V"}- We have 

d{ffn{u> -0}) = {Off n{u> ijj}) u n {n = ip}) , 

and 

u<g on S n {u > V}, 
u < ip on ^ D {u = ip} . 
By hypothesis ()3.4p . at least one of the two preceding boundary portions is non-empty. 

Case 1 (Boundary portion di^ \ T, {u > ip} is non-empty). Because u < g < v on dG \ S, then 
u<v on dG \ S n {n > ^/^} and we obtain n < u on ^ n {ti > ■;/'}• Thus, u <v on 6 . 

Case 2 (Boundary portion n {u = V} is non-empty). Because u = 'ip < v on G r\ {u = '\\)\^ 
then li < f on (9 (^ n {n > -0}) and we must have u<v on G r\\u> Thus, u < v on G. 

Hence, when G D {u > ^} is connected, we conclude that u < v on G. 

If G D {u > ip} is not connected, we may apply the preceding argument to each connected 
component C of ^ H {it > noting that G D {u > ip} is locally connectecfl and so each 
connected component of {n > V'}, in particular C, is a relatively open and closed subset. We 
have (see Remark l3.17p 

dC cd{Gn{u> ip}) = [pG r\{u> V'}) y^[Gr\{u = V'}) , 

and 

u<g on (9Cna^\Sn{u > ^}, 
on (9Cn {n = ^}. 

We now appeal to hypothesis (j3.4p and proceed exactly as before. This completes the proof. □ 

Remark 3.17 (The boundary of a connected component C of a subset X C M'^). This is a 
standard fact of point-set topology, but we include the proof for the sake of completeness. Recall 
that the maximal connected subsets (ordered by inclusion) of a nonempty subset X C M'^ are the 
connected components of X and are non-empty. Let C C X be a connected component and let 
X E (5, where C denotes the closure of C in M*^. Suppose that x E int(X). Since M'^ is locally 
connected, there is a connected set V d X which is an open neighborhood of x. Since C is 
connected and V r\C 7^ 0, then Cuy C X is connected. Therefore, V C C. That is, x E int(C). 
Thus, for any x E C, we have 

X E int(X) =^ x E int(C). 

Therefore, 

dC = C\ int(C) C C \ int(X) CX\ int{X) = dX. 
This proof works for any locally connected topological space in place of M"^. 

We can now state and prove the desired version of Proposition [331 where the requirement that 
the solutions or supersolutions belong to C'^{G) is removed. 

^For each point p £ C'r\{u > ip} and open subset V B p and V G ffr]{u> ?/)}, there is a connected open subset 
U withp G f/ C 1/ [55]. 
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Proposition 3.18 (Weak maximum principle and a priori estimates for solutions and superso- 
lutions to obstacle problems). Let G ^ and S ^ dff be as in Definition \3.6\ and let A in (jl.3p 
obey Hypothesis \3. 7\ and have the strong maximum principle property on ffUT,, and, in addition, 
require that c obeys (j2J3|) . Let f G C"(^), and g S Cioc(9^ \ S), and ip G Cioc(^) with ip < g 
on dG \ Yi. Suppose u is a supersolution to the obstacle problem in the sense of Definition \3.9\ 
and obeys (|3.4p . 

(1) lff>0 onG, then 

n > A inf q on G . 

(2) If f has arbitrary sign but there is a constant cq > such that c obeys (I2.14p . then 

n > A — inf / A inf 5f on G. 

Co di>\T, 

(3) If f < on G , and u is a solution for f and g and (Definition \3. 1\) . then 

n < V sup g V sup ip on G . 

(4) If f has arbitrary sign, u is a solution for f and g and ip, and c obeys ()2.14p . then 

u < V — sup / V sup g V sup ip on G . 

Co d0\T. 

(5) If u\ and U2 are solutions, respectively, for f\ > f2 and tpi > ip2 on G, and gi > g2 on 
dG \ Tj, and ui obeys ()3.4p . then 

ui > U2 on G. 

(6) If ui and U2 are solutions, respectively, for fi,ipi and /2,V'2 on G, and gi and g2 on 
dG\T,, and c obeys ()2.14p . andui,U2 obey ()3.4p . then 

hi - U2\\c{0) < — ll/i - h\\c{g') V \\gi -92\\c(d0^) V IIV'i - V'2||c(<f)> 

while if fi = / = /2 o.nd c just obeys (I2.13p . then 

\\UI - U2\\c(0) < Ibl - 92\\c(aff\^) V IIV'I - 1p2\\c{0)- 

Proof Consider Items ([I]) and When / < on i^, let 

m := A inf q, 

d0\T. 

while if / has arbitrary sign and A obeys (|2.14p . let 

m := A — inf f A inf q. 

Co d0\T. 

We may assume without loss of generality that m > — oo. Then m < g on dG \ S, while 

Am = cm < < / on G, 

when f < on G and 

Am = cm < CQm < f on G, 

when / has arbitrary sign. Hence, m is a subsolution and u a supersolution to the boundary 
value problem (II. ip in the sense of Remark 13.101 and m — li is bounded above on G (by Definition 
13. ip . Thus, m < n on 1^, as desired. 
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The proofs of Items ([3|) , ([3]) , ([5|) , and ([6]) here are the same as the proofs of Items ([3]) , (jH) , © , 
and dni) in Proposition 13.51 except that appeals to Proposition 13.31 are now replaced by appeals 
to Theorem 13.151 □ 

Remark 3.19 (A priori estimates for continuous supersolutions and subsolutions to a Dirichlet 
boundary value problem). Just as we extend Proposition 13.51 to Proposition 13.181 by allowing 
continuous rather than supersolutions, one could easily state and prove a version of Proposition 
12.191 where continuous supersolutions and subsolutions in the sense of Remark 13.101 replace 
supersolutions and subsolutions in the sense of Definition 12.141 



4. Strong maximum principle and applications to boundary value problems 

The usual statements of the Hopf boundary point lemma [401 Lemma 3.4] require that A in 
(|1.3p be uniformly elliptic on a domain, but a more careful analysis shows that it holds under 
much weaker hypotheses. We exploit our version of the Hopf lemma (see Lemma 14. ip to prove 
a strong maximum principle suitable for operators non-negative characteristic form (Theorem 
14. 6p and corresponding uniqueness results for solutions to equations with Neumann boundary 
conditions (Theorem 14.91 and Corollarv l4.12p . 

4.1. A generalization of the Hopf boundary point lemma to linear, second-order dif- 
ferential operators with non-negative characteristic form. Noting carefully that our sign 
conventions for differential operators and normal vectors are opposite to those of [30] and that 
we do not distinguish in the statement of lemma between points po S dff where a(po) = (that 
is, Po G S) or where a(po) is positive definite (that is, po G di^ \ ^)) we have 

Lemma 4.1 (Hopf boundary point lemma for linear, second-order differential operators with 
non- negative characteristic form). Suppose that <ZW^ is a domain and that G obeys an interior 
sphere condition at po G Off, with an open ball B{p,R) C ^ such that po £ dB{p,R). Require 
that the operator A in (11. 3|) obey (j2.ip on B(p,R), and 

c>0 onB{p,R), (4.1) 

and that at least one of the of the following conditions hold, where h = fi{po) denotes the inward- 
pointing unit normal vector at pq: 

(ah,h)>0 and > -2Ki and < Kq onB(p,R), (4.2) 

[ah, h) {ah, h) 

or OO and > i^o onB{p,R), (4.3) 

c 

for some non-negative constants, Kq,Ki, or some positive constant, vq, respectively. 
Suppose that u G C^(i^) obeys (j2.7p . that is, 

Au<Q on 0' 

and that u satisfies the conditions 

(i) u is continuous atpo; 

(ii) u{po) > u{p), for all p £ G; 

(iii) Dnu{po) exists, 

where Dnu{po) is the derivative of u at pQ in the direction of the inwaid-pointing unit normal 
vector, fi(po). Then the following hold. 
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(1) If c = on G , then Dnu{pQ) obeys the strict inequality, 

Dnu{po) < 0. (4.4) 

(2) Ifc>OonG and u{po) > 0, then (jM]) holds. 

(3) If u{po) = 0, then (|4.4p holds irrespective of the sign of c. 

Remark 4.2 (Comments on the Hopf lemma in the case of boundary points where the operator 
is degenerate). The Hopf lemma for points in {x2 = 0} for the Heston operator in Example 
11.21 was proved independently by P. Daskalopoulos in an unpublished manuscript using a barrier 
function similar to that in the proof of [40^ Lemma 3.4]. I am grateful to her for suggesting that 
such results should hold even at boundary points where the operator becomes degenerate. As 
pointed out to me by C. Pop, another version of the Hopf lemma for was obtained by Epstein and 
Mazzeo as \27\ Lemma 4.2.4] for their generalized Kimura diffusion operators, but also proved 
using a barrier function similar to that in the proof of [401 Lemma 3.4]. 

The traditional proof of the Hopf lemma, as described by Gilbarg and Trudinger \4:0\ Lemma 
3.4], by Evans [28\ §6.4.2], or by Han and Lin |44^ Theorem 2.5], exploits the interior sphere 
condition by choosing the barrier function 

^(^) = e-"l-'l' - e-°^', X G B{p,R) \ B{p,p), 

where a > is a constant which is ultimately depends on bounds on the coefficients of A on the 
open annulus B{p,R) \ B{p,p), where p € (0, i?) is a constant. This is the model for the barrier 
function chosen in the proof of \27\ Lemma 4.2.4], but the resulting argument is quite difficult. As 
we shall see in our proof of Lemma 14. H however, a choice of exponential- linear or linear barrier 
function instead will easily lead to the desired result. 

That such a Hopf lemma should hold even at points in S can be seen by examining the Kummer 
equation [H §13.1.1], 

Au{x) := —xuxx{x) — {b — x)ux{x) + cu{x) = 0, x £ M+, 

where b and c are positive constants here. If Ux{0) exists, then u is necessarily a constant multiple 
of a confluent hypergeometric function of the first kind, u{x) = kM{x), by [21 Sections 13.1.2-4, 
13.4.8, 13.4.21, and 13.5.5-10] and u is C°° on [0,00) with M(0) = 1 and u(0) = k£R. (Indeed, 
when a = c then M{x) = e^.) The continuity of Au on R_|_ implies that Ux{0) = kc/b. If A; < 0, 
then Au = on M+ and u has a strict local maximum at x = G M+ and Ux{0) < 0, as predicted 
by Lemma l4.ll 

Proof of Lemma \4-1\ We use the strategy of the proof of j40l. Lemma 3.4], with two choices of 
barrier functions, depending on whether condition (|4.2p or (14. 3p holds but both different from 
that in the proof of \40\ Lemma 3.4] and an approach to exploiting the interior sphere condition 
which is also different from that in the proof of [40 ^ Lemma 3.4]. 

Step 1 (Geometric set-up and application of a diffeomorphism) . We may assume without 
loss of generality, using a translation of M"^ if needed, that po = G M"^ and h = Cd and that 
B{p, R) is contained in the open upper half-space {xd > 0}. We now apply a diffeomorphism, 
^> : M'^ M'^ with $(0) = 0, to flatten the portion {0 < Xd < R} D dB{p, R) of the boundary of 
B{p, R) by pushing it downward, so that 

^ ({0 <Xd<R}n dB{p, R)) = {x£R'^ : \x\ <R,Xd = 0} 

= BiO,R)n{xd = 0}. 
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If T C di? is a small, relatively open neighborhood of po £ the map <I> pushes T \ {po} 
downward into the open /ower half-space, {x^ < 0}: 

$(r \ {po}) c {xd < 0} 

> R}ndB{p,R)) ^ $(^). 

Henceforth, after applying the preceding diffeomorphism (and now denoting $(^) simply by ff), 
we may assume, without loss of generality, that the open half-ball 

B+{0, R) := {xeR'^ : \x\ <R,Xd>0} C 

has the property that 

{xd > 0}ndB+{0,R) (E ff. 
One can easily check (see the proofs of [281 Theorem 6.3.4], |401 Lemma 6.5 or Theorem 8.12], 
\4:8\ Lemma 6.2.1] for similar arguments) that the conditions (14. 2p and (]4.3p . thanks to (j2.ip . 
are invariant under the preceding coordinate changes, with the only modification being that the 
non- negative constants Kq,Ki and positive constant i^q may be replaced by different non- negative 
constants or a positive constant, respectively, and the open ball, B{p,R), in those conditions is 
replaced by the open half-baU, B^{0,R). 

Step 2 (Construction of the barrier function when condition (14. 2p holds). We choose 

v{x) := e"^'* - 1, xe 

where a > is a constant yet to be determined. Clearly, v{0) = and v > on the half-space 
{xd ^ 0} and, in particular, v > on B^{0,R). Moreover, 

Av = -a^a'^'^e''''^ - + c (e°^'* - 1) 

< -a'^'^ (a^ + 2Kia- Kq) e"^'* on S+(0, R), 

using the hypothesis (|4.ip to obtain c > on B~^{0, R) and noting that a'^'^ > on -6^(0, R) by 
hypothesis (|4.2p and that a,b,c obey (j4.2p (after applying a diffeomorphism) with B{p,R) 
replaced by B+{Q,R). But 

+ 2Kia - Kq = {a - Kif - Kf - Kq > 

provided a obeys 



a>Ki + ^Ko + Kf. 
We fix such an a and thus obtain Av < on B^(0, R). 

Step 3 (Construction of the barrier function when condition (14. Sp holds). We choose 

v{x) := Xd, Vx G M'^, 

and observe that 

Av = -h'^Vx^ +CV = - f^ + cxd- 
We may assume without loss of generality, by decreasing R if needed, that B~^{0, i?) C {x £ M*^ : 
< Xrf < vq} for fo as in (|4.3p . and so the bound (j4.3p (after applying a diffeomorphism) with 
5(p,i?) replaced by B+{0,R) yields 

inf — > 1^0 > 0. 

B+{0,R) C 
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Therefore, 

Av = -cl Xdj < -c(fo - Xd) < on B~^{0, R), 

and we again obtain Av < on -6^(0, R). 

Step 4 (Verification tliat tlie weak maximum principle holds for A on B'^{0, R)). From Steps [2] 
andO we obtain Av < on -B+(0, R) for either choice of barrier function, v. Therefore, because A 
obeys (j2.ip on B~^{0, R) (note that this condition is invariant under the changes of coordinates 
employed here), the proof of |401 Theorem 3.1], with the role of e'''^^ in [lOl p. 32] replaced by v, 
shows that the conclusions of [40 1, Theorem 3.1 &: Corollary 3.2] hold for A on B^{0,R). 



Step 5 (Application of the weak maximum principle). Since u — u{0) < on ^ and u G C\o 

and {xd > 0} n dB+{0, R) d ff, we obtain 

u{x) - n(0) < -mo < 0, Vx € {xd > 0} n dB+{0, R), 
for some positive constant, ttiq, depending on R and u. If (14. 2p holds, then 

v{x) = e"^'* - 1 < e"^ - 1, Vx G {xa > 0} n dB+{0, R), 
while if ()4.3p holds, then 

v{x) = Xd<R, yx £ {xd> 0}ndB+{0,R). 
Hence, recalling that R > 0, there is a positive constant mi := (e°^ — 1) V i? such that 

v{x)<mi, yx e {xd> 0}ndB~^{0,R). 

Consequently, 

u{x) - n(0) + ev{x) < -mo + emi < 0, Vx G {x^ > 0} n dB~^{0, R), 

provided we fix e in the range < e < mo/mi, while 

n(x) - n(0) + et;(x) = u{x) - u{0) < 0, Vx G {xd = 0} n 95+ (0, i?), 

since for either choice of barrier function we have v{x) = when x^ = and our hypothesis (jn]) 
(with po = 0) implies that 'u(x) < u{0) on dB+{0,R) C {0}. But 

A{u - u{0) + ev) =Au- cu{0) + eAv < -cn(0) < on 5+(0, R), 

where the last inequality holds if c = (as in Conclusion ([T])), or c > and u{0) > (as in 
Conclusion ([2])), or u{0) = (as in Conclusion ([3])). (For the case n(0) = 0, we simply note as in 
the proof of [40^ Lemma 3.4] that we can replace Ahy A + c~ , where we write c = c"*" — c~.) 
The weak maximum principle therefore yields 

u-u{0)+ev<0 onB+(0,i?), 

by virtue of Step HI 

Step 6 (Sign of the directional derivative of the subsolution at the boundary). We have 
^iM^ < ^ VxGi?+(0,i?), 

Xd Xd Xd 

If v{x) = e"^** — 1, we have = ae"^'' and fxjj(O) = a > 0, while if v{x) = Xd, we have Vx^^ = 1. 
Taking the limit 4 and noting that 

a if g. b. c obey (|4.2p . 
1 if 6, c obey ([13]), 



PARTIAL DIFFERENTIAL OPERATORS WITH NON-NEGATIVE CHARACTERISTIC FORM 



29 



yields Vx^ (0) > a A 1 and 

and thus (jili]) holds. 

This completes the proof. 



UxM < -eVxM < -eia A 1) < 0, 



□ 



Remark 4.3 (Application to the elliptic Heston operator). The hypotheses of Lemma 14.11 on 
the coefficients of A are obeyed in the case of the elliptic Heston operator, Example 11.21 where 
d = 2 and 0' G M and S = mt{dff H dM). For example, po £ Ti then h = n(po) = ^2, while 



,22 



and = k{6 - y), so 



2K{e - y) 



dP- a^y 

Thus, condition (j4.2p is obeyed when r = 0, noting that c = r, while if r > 0, then 

62 K{e - y) kO 



> 



2k , , 
— , Vy>0. 



>—, 0<y<e/2, 
2r 



and thus condition ()4.3p is obeyed. 



□ 



Remark 4.4 (Application to the linear, second-order, uniformly elliptic operators). Given a 
domain C M"^, if A in ([O]) obeys (f2T]l . one has 



< A(a;) < {a{x)r], r]) < A{x), Vx G ^, r? G M'^ \ {0}, 

where X{x) and A(x) are the minimum and maximum eigenvalues of a{x). Then, A is elliptic on 
0' in the sense of [401 p. 31] if A > on ^, and uniformly elliptic on 0' |40l p. 31] if 



sup 



< oo. 



while A is strictly elliptic on j4Ul p. 31] if 



inf A > 0. 



(4.5) 



(4.6) 



(The terminology is not universal — for example, compare |28t §6.1.1] — but we shall follow the 
convention of |40j in this article.) The hypotheses of [40^ Lemma 3.4] require that A is uniformly 
elliptic on ^ and that (see [iQl Equation (3.2)]) 



sup 



< oo and sup 

0- 



< OO. 



(4.7) 



The bounds in (j4.7|) . together with an assumption that A is elliptic on ff, ensure that the in- 
equalities ()4.2p hold and hence our Lemma l4. 1 1 implies |40^ Lemma 3.4], although the converse is 
not true, as illustrated in Remark 14.31 □ 



Remark 4.5 (Hopf boundary point lemma for domains obeying an interior cone condition). The 
interior sphere condition can be relaxed in the classical Hopf lemma \4:0\ Lemma 3.4], as noted 
in \40\ p. 35 & p. 46], and generalizations to domains with non-smooth points are described in 
\52 \ \53 \ [56 l 157]. We shall describe a generalization of Lemma l4.1l to such domains in |29j . □ 
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4.2. Strong maximum principle. By using Lemma im instead of [iQl Lemma 3.4] in the proofs 
of |40| Theorem 3.5 & 3.6], we obtain the following statement of the strong maximum principle 
in Theorem 14.61 We provide details of the proof in order to emphasize the fact that the phrase 
"maximum in the interior of i^" in our version of the strong maximum principle can be replaced 
by "maximum in S or the interior of following the pattern observed by Daskalopoulos and 
Hamilton in |19j where S can often be treated as a set of "interior" points of the domain. 

Theorem 4.6 (Strong maximum principle for functions). Suppose that G C M"^ is a domai?^. 
Require that the operator A in ()1.3p obey (j2.ip . ()2.3p . ()2.13p . and that on each ball B C 0' with 
Be U either ()4.2p or ()4.3p hold, where S C di^ is as in ()2.2p . Assume, in addition, that 
S is and 

{b,n)>0 onT,. (4.8) 
Suppose that uGC'^{ff)n Cioc(^) obeys (1231), and ([Sj]). If u obeys ^J^, that is, 

Au<Q on 

then the following hold. 

(1) If c = on 6 and u attains a maximum in S, then u is constant on G . 

(2) If c> on G and u attains a non-negative maximum in S, then u is constant on G . 

Remark 4.7 (Strong maximum principle property). Theorem 14.61 implies that A has the strong 
maximum principle property, for functions, in the sense of Definition 13.141 

Proof of Theorem \4-6[ Consider Conclusion ([T]). Assume, to the contrary, that u is non-constant 
on G and achieves its maximum M on at a point in i^U S. Let := {x £ G : u{x) < M} 
and observe that is non-empty by our assumption that u is non-constant on G. Let p S 
be such that dist{p,dff~) < dist{p,dG \ S) (if S = dff, then any p e G~ wih do) and let 
B C be the largest open ball centered at p and contained in . Then u{po) = M for some 
PoedBn Off- and u < M on B. Note that po^ ffUT., since dist(p, d^') < dist(p, \ S) by 
construction (again, this condition is vacuous if S = dff). 

Case 1 (pq £ G). We must have Du{p()) = since pQ is an interior local maximum. However, by 
applying Conclusion ([T]) in Lemma 14.11 to the operator A on the domain B and boundary point 
Po G dB, we obtain Du{po) 7^ 0, a contradiction. 

Case 2 (po G 5]). We may assume without loss of generality, by applying a diffeomorphism 
if needed, that ^ n B{pQ,r) = H n B{pQ,r) for sufficiently small r, where IHI = M'^"^ x M+, so 
S n B{po,r) = {xd = 0} n B{po,r) and n(po) = e^. Using ^M), and b'^{po) > by (jilH]) . and 
Uxiipo) = for 1 < i < d — 1 since po is a local maximum for u in E, and the fact that c{po) > 
and u{po) > by hypothesis, we would then have 

Au{po) = -tT{aD'^u){po) - {b{po),Du{po)) + c{po)u{po) > -b'^{po)uxa{Po) > 0, 

where the final strict inequality follows from Conclusion ([TJ of Lemma 14.11 contradicting the fact 
that Au < on ^ by hypothesis and hence Au < on ^ U S by (j2.5p . 

Conclusion ([2]) follows by an identical argument when u achieves a non-negative maximum in 
i^U S, except that we now appeal to Conclusion ([2]) in Lemma |4. 11 □ 

Theorem 14.61 provides an interesting conclusion when S = di^ and is bounded: 



Recall that by a "domain" in R'', we always mean a connected, open subset. 
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Corollary 4.8 (Strong maximum principle for functions on bounded domains and operators 
which are degenerate along the entire domain boundary). Assume the hypotheses of Theorem \4-6\ 
and, in addition, that T, = dG and G is bounded and c = 0. Then u is constant on G . 

Proof. Since u G C{G) and u is bounded, then u necessarily achieves its maximum at some point 
in ^ = ^US. But then Conclusion ([T|) in Theorem 14. 61 implies that u must be constant on G. □ 

We next consider the question of uniqueness in the Neumann problem and note here that it is 
important to distinguish between dG and dG \ S. 

Theorem 4.9 (Uniqueness for the Neumann problem for functions on bounded domains). 
Let G C be a bounded domain. Let A in (jl.3p and S ^ dG in (j2.2p obey the hypotheses 
of Theorem \4-6\ and assume that G satisfies an interior sphere condition at each point of the 
boundary portion dG \ S. Suppose that u G C'^{G) n C{G) obeys (pli]) . (p3]) . ([Z^G]) and that 

An = on G. 

If the derivative, DnU, with respect to the inward-pointing normal vector field, n on dG \T,, is 
defined everywhere on dG \ S and 

Dnu = ondG\T., (4.9) 

then u is identically constant on G. In addition, if c > at some point in G, then n = on G. 

Proof. If u is not identically constant on G, then either u or —u achieves a non- negative maximum 
M on G. Since G is bounded and u £ C{G), we may suppose that u achieves a non-negative 
maximum at some point pq G G, as the argument when —u achieves a non-negative maximum on 
G will be identical. Therefore, u{pq) = M for some po S dG\T, since, because u is not identically 
constant on G, Theorem 14.61 implies that u < M on G L)Y1. But then Lemma l4.ll implies that 
Dnu{po) < 0, contradicting our hypothesis ()4.9p . Thus, we must have u = M, a constant, on G. 
If c > at some point of G, the facts that Au = cM and Au = force M = 0. □ 

Remark 4.10 (Interior sphere condition in the hypotheses of Theorem 14. 9p . In applications to 
mathematical finance, the hypothesis that G obey an interior sphere condition at each point 
of the boundary portion dG \ S is not usually obeyed at points in dT, C dG \ S. Alternative 
hypotheses, more in harmony with applications, are discussed in |29j . 

Definition 4.11 (Classical solution to a boundary value problem with Neumann data). Let 
G C M"' be a domain, let S ^ dG, and require that G satisfy an interior sphere condition at each 
point of the boundary portion dG \ S. Given functions / G C{G) and h G C\oc{dG \ S), we call 
u G C'^{G) n C\oc{G) a solution to a boundary value problem for a linear, second-order partial 
differential operator A in (jl.3p with Neumann boundary condition along dG \ S, if u obeys ()2.4p . 
()2.5p . ()2.6p . and the derivative, DnU, of u on dG \ S in the direction of the inward-pointing unit 
normal vector field, n, is defined everywhere on dG \ S, and 

Au = f on G, (4.10) 

DnU = h on G\T.. (4.11) 

Corollary 4.12 (Uniqueness for the Neumann problem for functions on bounded domains). 
Let G C R'^, and A, and f e C{G), and h e Cioc{dG \ S) be as in Definition \4-Il\ Require in 
addition that G be bounded and that c > at some point of G\Ti. If ui,U2 are solutions in the 
sense of Definition \4-Il\ then ui = U2 on G. 

Proof. By hypothesis, A{ui — U2) = on ^ and Dn{ui — U2) = on dG \ S, so ui — ii2 = on 
^ by Theorem 1121 □ 
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5. Weak maximum principle for smooth functions 



Having considered applications of the weak maximum principle property (Definition I2.8P to 
Dirichlet boundary value problems in fJ21 obstacle problems in fJS]and, in conjunction with a Hopf 
Lemma for operators with non-negative characteristic form, to the strong maximum principle and 
Neumann boundary value problems in 21 we now establish conditions under which the operator 
A in (ll.3p has the weak maximum principle property on U S. In §5.H we establish a weak 
maximum principle for bounded functions on bounded domains (Theorem 15. ip . while in Hb.2\ 
we extend this result to the case of bounded C'^ functions on unbounded domains (Theorem 15. Sp . 

Our weak maximum principle (Theorems 15.11 and 15 . 3p differs in several aspects from [631 Theo- 
rem 1.1.2], some of which may appear subtle at first glance but which are nonetheless important 
for applications. For example, 

(1) The function u is not required to be in C^{ffUT.) fl C(^), but rather C'^{ff) n C(^) and 
obey (EaD, (I23|), (USD; 

(2) The subdomain C M'^ is allowed to be unbounded; and 

(3) The coefficient functions of the partial differential operator A in (|1.3p are allowed to be 



The significance of these points is illustrated further by the example of the Heston operator 
discussed in Appendix ICl 

5.1. Bounded functions on bounded domains. We begin with the case of bounded 
domains and adapt the proofs of [401 Theorem 3.1] and |481 Theorem 2.9.1]; see also |191 Theorem 
1.3.1], |30| §3.2]. It will be convenient to adopt the following convention. If S ^ di? and 
g : dG \ S — )• M is a function and m E M, then 



where we recall that x V y = max{x, y}, for any x, y € M. 

Theorem 5.1 (Weak maximum principle for functions on bounded domains). Suppose that 
G <Z is a bounded domain, with S as in Definition \2.1\ Require that the coefficients of the 
operator A in (jl.3p be defined everywhere on G , obey ()2.ip an(0 



unbounded. 




(5.1) 




where S is assume^^ to be with inward-pointing unit normal vector field n. Assume further 
that at least one of the following holds, 



c> on S, or 



S ^ dG and {b, h) > on T, and inf 




> — oo 



(5.4) 



for some fixed direction, h G M'^, with \h\ = 1. 



The condition on c here is shghtly stronger than that of (|2.13|) . 

'This assumption is required by the change-of-coordinates argument used in Step [1] of the proof. 
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Suppose that u G C'^{G)r\C{G) obeys (H^j), and If u obeys ([22]), that is, ylu < 

on G , then 

supii < V sup u, (5-5) 



and, if c = on GUT, ( and so ^ dG by our hypothesis (j5.4p ), then 



supu = sup u. (5-6) 



Moreover, A has the weak maximum principle property on GUT, in the sense of Definition \2. 

Remark 5.2 (Application to the elliptic Heston operator). The hypotheses of Theorem 15.11 on 
the coefficients of A are obeyed in the case of the elliptic Heston operator, Example 11.21 where 
d = 2 and ^ C IH is bounded and E = ^ n dM.. Choosing /i = 62 in condition ()5.4p . we have 
6^(x) = k{0 — X2) and a'^'^{x) = 0-^x2/2, so 

(6,62) b'^ix) 2K{e-X2) 2^6 2k 2k ^ 

— — — —^> ^) V3;2 > 0, 



and condition ()5.2p is obeyed (with n = h) as (6, /i) = ^^(xi, 0) = k0 > 0. Lastly, c = r > 0. □ 

Proof of Theorem \5.1[ Since u is continuous on G and by hypothesis ^ is compact, then u achieves 
its maximum value at some point p in G. 

Step 1 (c > on GUT,). Suppose that u attains its maximum value, u{p), at a point p in GUT,. 
If p e G, then Du{p) = and D'^u{p) < 0, so that (j2T]) gives 

An(p) = -a'^ {p)ux,xj{p) - b\p)ux^{p) + c{p)u{p) 

> c{p)u{p), 

since a(p) is non-negative definite and tr{KL) > whenever K, L are non-negative definite 
matrices [50l p. 218]. If u{p) > 0, we would obtain Au{p) > 0, contradicting our assumption 
()2.7p that Au < on G. Therefore, we must have u{p) < and, necessarily, u{p) < supg^yj] n"*" 
when dG \ S is non-empty. 

If p G S then, possibly after a change of coordinates on M*^, we may assume without loss 
of generality that B{p) n{GUJ:) = B{p) n H, where H = M'^"^ x M+. Thus, n{p) = ea and the 
condition (j5.2p at p becomes b'^{p) > 0. We have tT{aD'^u{p)) = by ()2.6p . while Ux^ip) = for 
1 < i < d — 1 and Ux^{p) < since p G S is a local maximum and ()2.4p holds, so that 

Au(p) = -a'-^{p)ux,xj{p) - V{p)uxi{p) + c{p)u{p) 
= -b'^{p)ux^{p) + c{p)u{p) 

> c{p)u{p). 

Our assumption (12. 7p that Au < on G implies, by continuity of Au on GuT, via (|2.4p and (|2.5p . 
that Au < on i^U S. If u{p) > 0, we would obtain Au{p) > 0, a contradiction. Therefore, we 
must again have u{p) < and, necessarily, u{p) < supg^'^j] u"^ when dG \ S is non-empty. 

Finally, if p G dG \ S, then u{p) = supg^^^ u and so by combining the preceding three cases 
we obtain (jS.Sp for this step. 
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Step 2 (c = on ^ U S an(£3 H ^ diff and Au < on ^ U S). Suppose that u attains its 
maximum value, u{p), at a point p in S. Repeating the argument of Step [1] would then yield 

Au{p) > c{p)u{p) = 0, 

contradicting the assumption that Au < on i^. Consequently, we must have p E d& \ S and 

supii = sup 

that is, (|5.6|) holds. 

Step 3 (c = on U S and S ^ dff). We see that, for a constant > to be chosen and any 



\ a^^ [x)hihj J 



b^(x)hi 

< on ^, if inf ..; /, \ > -ly, 
x£ff a^^[x)hinj 

appealing to (j5.4p for the ability to choose a finite v > For any x G S, we have 

Ae'^'^^''^) = -iyb\x)hie''^'''''^ 

< on S, if b'hi > on S, 
again appealing to (j5.4p . Therefore, we have, for any e > 0, 

Aiu + ee"^'"'''^) <0 on^US, 

and so, by Step [H we obtain 

sup(u + ee^<'^'^>) = sup {u + ee'^<'^'^>). 
Letting e — )• yields (j5.6p for this step. 



Step 4 (c > on ^ U E and S ^ dff). Let ff^ denote the open subset {p e : u{p) > 0} C 
If is empty, then u < on ^ and so u < on ^, with 

sup n < = sup . 



It remains to consider the case where is non-empty. By hypothesis (12. 7p . j4ti < on and 
thus Aqu < —cu < on where Aq := A — c. We may write 

= (s n u \ S) . 

If (9(^+ \ S were empty, then we would have dff~^ C S C c?i^. Thus, we would necessarily have 
0"^ = G and hence = S = dG ^ contradicting our assumption that S ^ dG . Therefore, 

\ S must be non-empty and our result (15. 6p for the case c = on U S yields 

sup u = sup u. 



^^If E = d6, then the conditions < on E and tr(aD^u) = on E via (|2.6|1 imply that -(b, _Du) < on E. 
If in addition we had h = {h,n)n and (|5.2|l were strengthened to (b, n) > on E, then we would obtain D„u > 
on E and consequently u would have a local maximum in ff, contradicting the condition Au <Q on ff. In general, 
when E = dff, one obtains no additional information regarding sup^ u when c = on ^ U E, beyond the fact that 
u achieves its maximum at some point of U E = ^. 
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Now 

\ s = n a€? \ s) u n ^) . 

Since u = on 9^"^ n & (because u is continuous on ^ and ti < on \ G^)^ then 

sup n = V sup u = sup n, 

9<?+\S d0+r\ae\T. d0\Y, 

and so, combining the preceding inequahties and equations and noting that sup^ u = sup^+ u, 
we obtain 

sup u = sup u, 

as desired for the case of non-empty Combining the preceding cases yields ()5.5p for this step. 
This completes the proof. □ 

5.2. Bounded functions on unbounded domains. Next, we consider the case of bounded 
functions on unbounded domains. We have the following refinement of the maximum prin- 
ciple for bounded functions on unbounded domains and elliptic operators with non-negative 
characteristic form [48 1. Theorem 2.9.2 &: Exercises 2.9.4, 2.9.5]. 

Theorem 5.3 (Weak maximum principle for bounded functions on unbounded domains). 
Assume that the coefficients of A in (11. 3p obey the hypotheses of Theorem \ 5.1\ except that the 
condition (|5.4p on (a, b) or c is replaced by the condition (j2.14p on c for some positive constant, 
Co, and, in addition, we require that there is a positive constant, K , such that 

tra(x) + (6(x),x) < K(l + Ixp), Vx G ^. (5.7) 

Suppose that 0' ^M.'^ is a possibly unbounded domain and that u G C^(^) H C\oc{^) obeys (12. 4p . 
(12. 5p . and (|2.6p . //sup^n < oo and u < on d& \ S (when non-empty), then 

1 

sup u < V — sup Au. 

e Co 

Moreover, A has the weak maximum principle property on £?U S m the sense of Definition \2.8[ . 
Proof. Let 

t;o(x) := 1 + VxGM'^, (5.8) 

and observe, noting that c > on ^ by hypothesis (j2.14p . that 

Avo{x) = -2tra(x) - 2{b{x),x) + c(x) (l + |xp) 
> -2Kvo{x), Vx G ff, 

where K is the constant in ()5.7p . Therefore, 

(A + 2K)vo > on ^. (5.9) 

Define 

M := V sup(^ + 2K)u. 

e 

If M = +00, there is nothing to prove, so we may assume < M < cxo. For 5 > 0, set 

u; := ti - (5?;o - (co + 2A0~^M, (5.10) 
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and observe that, using c > cq > 0, 

{A + 2K)w = {A + 2K)u - 5{A + 2K)vq - (c + 2K){co + 2K)-^M 
<{A + 2K)u - M 
< on ^. 

Denoting B{R) = {x S : \x\ < R}, then for ah large enough i? > 0, we have < on 
^ n dB{R) since u is bounded above on ^ by hypothesis. Also, u < on 9^ \ S (when non- 
empty) by hypothesis and so w < on 9^ \ S (when non-empty). But 

d{ff n B(i?)) = (^ n as(fl)) u {B{R) n a^), 

and therefore, 

d{ff n B(/?)) \ s = (^ n dB{R)) u (S(i?) n \ s), 

so that 

w<0 on 9(^nS(ii)) \S. 

Consequently, Theorem l5.ll (with A replaced by A + 2K) implies that w < on 0'nB{R), for any 
sufficiently large R > and thus < on ^. By letting (5 — )■ 0, we obtain u<{co + 2Ky^M on 
0" and so 

, M 
supti < 



and thus 



C0 + 2K 
= sup((A + 2K)n)+ 

< sup + . supu^ 

Co 2A CQ + 2K 



Co I 1 

■supn < -— sup(^ti) 



Co + 2K / - Co + 2K 
which, using cq > 0, completes the proof. □ 

Remark 5.4 (Application to the elliptic Heston operator). The hypotheses ()2.14p and (j5.7p in 
Theorem 15.31 are obeyed in the case of the elliptic Heston operator (Example II. 2|) with 

aix) = Y (^pfj ' = {r-q- X2/2 Hi{e - X2)) , and c{x) = r, 

provided r > 0. □ 



Part 2. Weak maximum principles for bilinear maps and operators on functions in 
Sobolev spaces and applications to variational equations and inequalities 

In this part of our article (^|6] and ^J7|), we develop weak maximum principles for bilinear 
maps and operators on functions in Sobolev spaces and applications to variational equations and 
inequalities. 
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6. Applications of the weak maximum principle property to variational equations 

Just as in the case of the weak maximum principle for linear, second-order, partial differential 
operators A in (|1.3p with non-negative characteristic form acting on smooth functions, we shall 
encounter many different situations (bounded or unbounded domains C M'^, bounded or un- 
bounded functions u with prescribed growth, and so on) where the basic weak maximum principle 
holds for bilinear maps o on H^{0',K)) or associated linear operators A : iJ^(^, tt)) — )• L'^{ff,n). 
Again, we find it useful to isolate that key property and then derive the consequences which 
necessarily follow in an essentially formal manner. In this section, we consider applications to 
variational equations. After providing some technical preliminaries, we proceed to the main appli- 
cations, including the comparison principle (Proposition 16. 7( ) and a priori estimates (Proposition 
I6.10p for i^^((^, tr) supersolutions and solutions to variational equations, and the corresponding 
results (Proposition 16.15]) for ff^(^, tt)) supersolutions and solutions to the associated boundary 
value problems. In §2.31 we show that when a bilinear map o on H^{ff,to) or operator A on 
H'^{ff,to) has a weak maximum principle property for subsolutions (on unbounded domains) 
which are bounded above, the property may extend to subsolutions which instead obey a growth 
condition (Theorem 16.161 and Corollary 16. 17p . 

Definition 6.1 (Weight function). Let ^ C Rrf be a domain. We call tt) a weight function if 

W G C7(^) n L^^) and tt) > on ^. (6.1) 

Definition 6.2 (Weighted Sobolev spaces). Given weight functions tt>k,i, for integers < i < k, 
we define Hilbert spaces with norms 

IMH^iff^rv) X] / \D'u\'^n3k,idx, k>l, 

as the completions of the vector space C^{^) with respect to the preceding norms and denote 
L2(^, tt)) := H^{ff, tt)). Given S g d^, we define H^{ff U S, tt)) to be the closure of C^{ff U S) 

When S = 5^ in Definition ESI then ^ U S = ^ and C^{ff U S) = C^iS") and H^{^, tt)) = 
H^{ff,m). By analogy with [H §8.1], for u E H\ff,Ki), we define 

ess supu := ini {I £ R : u < I on dff \ T, in the sense of i7^(^, tt))| , (6.2) 

where we recall that u < I on dff \ T, in the sense of H^{0', tt)) if (n - /)+ G H^{0' U S, tt)). If 
g £ H^{ff, to) and m G M, then 



ess sup g if S ^ 

m V ess sup g ■= I dff\j: (6.3) 

m if s = a^. 



by analogy with our convention (j5.ip for everywhere-defined functions. 

Definition 6.3 (Weak maximum principle property for a bilinear map). Let ^ ^ M'^ be a domain, 
let 
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be a bilinear map, and let S ^ dG be a relatively open subseiF^. We say that o obeys the weak 
maximum principle property on U S if whenever u G //^(^, ro) obeys 

a{u,v) <0,\fv e H^{GL}T,,ni) with v >0 a.e. on ff, 
M < on 9^ \ S in the sense of H^{ff, tt)), 

and, if ^ is unbounded, ess sup^ u < oo, then 

u < a.e. on i^. 

Remark 6.4 (Auxiliary boundedness condition in the definition of the weak maximum principle 
property). Once a is known to have the weak maximum principle property on ^US, the condition 
that u be essentially bounded when ^ is unbounded (which could be replaced by a growth 
condition) plays no role in most applications of the weak maximum principle property. □ 

Remark 6.5 (Positivity conditions required for a bilinear map to have the weak maximum 
principle property on ^US). We do not explicitly require that a(l,-) > (note that 1 e iJ^(^,tt)) 
by (j6.ip ). but we impose this requirement when needed in addition to our assumption that o has 
the weak maximum principle property on i^US (compare |40( Equation (8.8)]). In order to prove 
that a has the weak maximum principle property on U S when is unbounded, we shall need 
to assume in addition that o obey the positivity condition (16. 9p . □ 

Examples of bilinear maps, a, with the weak maximum principle property on ^US are provided 
by Theorems lO K8[ KTR KTR [6461 and (when T, = 0) ^ Theorem 8.1]. 
By analogy with j28l §5.9.1], we let 

i7"^(^US,tt)) := (//(J(^UE,tr))' (6.4) 

denote the dual space and, as in [281 §5-9.1], observe that 

i/oH^US,ra) C L^{G,to) C i7-^(^UE,tt)). 

The proofs of [3l Theorem 3.8] or |28[ Theorem 5.9.1] easily adapt to show that every F € 
H~^{ffU T.,n) has the form 

F{v) = {f,v)L2(^^^^^-{f\v,jL2(^^^^), Vt;G Fi(^US,tt)), (6.5) 

and we write F = {f, f^,. . . , f^), where /, /* G L'^{ff, tt)), 1 < i < d. We say that F < when 

F{v) < 0, G H^{GUT.,to),v > a.e. on ff, 

and, given Fi, F2 e H-^{G U S, ro), we say that Fi < F2 if Fi - F2 < 0. 

Definition 6.6 (Variational solution, subsolution, and supersolution). Let F = (/, . . . , f^) € 
H~^{i^U S,tt)) and g G {iff, to). We define u G iJ^(^, tt)) to be a variational subsolution, 

a{u, •) < F, 

u < g on 80' \ S in the sense 

if ess sup^i u < 00 when 0' is unbounded and 

a{u,v)<F{v), yv € H^{0UT,,xo),v>O a.e. on ff, (6.6) 

(n-g)+ G//o'(^US,tt)). (6.7) 

^^We shall assume that E ^ dff is relatively open for the sake of consistency with Part [T] of our article, although 
E could now be any non-empty subset, not necessarily relatively open. However, the assumption of relative openness 
involves no loss of generality. 
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We call u E H^{i^,ni) a variational supersolution if —u is a variational subsolution and call 
u E H^{ff, vo) a variational solution if it is both a variational subsolution and supersolution. □ 

The first application of the weak maximum principle property, as in the case of Proposition 
12.161 is to settle the question of uniqueness. 

Proposition 6.7 (Comparison principle for variational supersolutions and uniqueness for vari- 
ational solutions). Let F E H^^{ffU Ti,tv), and g E H^{ff,tv), and a be a bilinear map on 
H^{ff,n) obeying the weak maximum principle property on ^ U S, for some T, ^ dff. Suppose 
that u E H^{i?,n)) is a subsolution and U2 E H^{ff,tt)) a supersolution in the sense of Definition 
\6.6l Then U2 > ui a.e. on G and ifui,U2 are variational solutions, then U2 = ui a.e. on G . 

Proof. Since u\ is a subsolution and U2 a supersolution, then u\ — U2 is a subsolution with 
ui — ii2 = on dG \ Yj in the sense of H'^{G, id) and thus 

lii — U2 < a.e. on G, 

because a has weak maximum principle property on S. When ui,U2 are both solutions, then 
we also obtain U2 — ui < a.e. on G and U2 = ui a.e. on G. □ 

We shall occasionally need the following analogue of Hypotheses 12.171 and 12.181 respectively. 

Hypothesis 6.8 (Non-negative zeroth-order coefficient). If o is a bilinear map on H^{G',K)), 
then 

a{l,v)>0, E i?o(^US,tt)),v > a.e. on ^. (6.8) 

Hypothesis 6.9 (Uniform positive lower bound on zeroth-order coefficient). If a is a bilinear 
map on H^^G, tt)), then there is a constant cq > such that 

a{l,v) > {co,v)]^2(^^'^n), \/v e H^{GuJ:,to),v >0 a.e. on G. (6.9) 
We can now proceed to give the expected a priori estimates. 

Proposition 6.10 (Weak maximum principle and a priori estimates for H^{G, ro) functions). Let 

F = (/, /I, . . . , /^) G H-^{& U S, id), and g E H^{&, Vo), and a be a bilinear map on H^{ff, ro) 
obeying the weak maximum principle property on ^US, for some S ^ di?, and assume that (16. 8p 
holds. Suppose that u E if^(^, ro). 

(1) If F <0 and u is a variational subsolution for F and g ( Definition 1 6. 6\) . then 

n < V ess sup 5 on G . 

(2) If F = {f,0, . . . , 0) and f has arbitrary sign and u is a variational subsolution for F and 
g but, in addition, there is a constant cq > such that a obeys ()6.9p . then 

u < V — ess sup / V ess sup g on G . 

Co d0\T. 

(3) If F >0 and u is a variational supersolution for F and g (Definition \6. 6)) . then 

u > A ess inf q on G . 

(4) If F = {f,0, . . . , 0) and f has arbitrary sign, u is a variational supersolution for F and 
g, and a obeys (|6.9p . then 

u > A — ess inf f A ess inf q on G . 

Co G de\T. 
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(5) If F = and u is a variational solution for F and g (Definition \6. 6)) . then 

\\u\\L°^{ff) < \\9\\L°°[d&\T,)- 

(6) If F = {f,0, . . . , 0) and f has arbitrary sign, u is a variational solution for F and g, and 
a obeys (|6.9p . then 

< — V ||fl'llL°°(9<5'\S)) 

The terms ess supg^x^-^g, and ess infg^^g 5; ^.''^d \\g\\L°°[d^\T,) the preceding items are omitted 
when T, = dff. 

Proof. The proof follows almost the same pattern as that of Proposition 12.191 For Items ([TJ and 
(l2|), we describe the proof when S ^ d^; the proof for the case S = dff is the same except that 
an essential supremum of a non-negative function over 80' \ S is replaced by zero. When / has 
arbitrary sign, choose 

M := V — sup / V ess sups', 

CQ f? diff\T. 

while if / < a.e. on 0, choose 

M := V ess supg. 

We may assume without loss of generality that M < oo. Since M > on dff \ S (in the sense 
of tt))) and cqM > / a.e. on ff, then, for all u G U S, ro) with u > a.e. on ff, the 

condition (|6.9p gives 

a{M,v) > icoM,v)L2(^ff^^) > (/,t;)L2(^,„), yv G H^{0U^,m),v > a.e. on ff, 

when / has arbitrary sign and, when / < 0, the condition (16. Sp gives 

a{M,v) > > {f,v)L2(^ff^ro), G i7o^(^US,rD),u > a.e. on 0. 

Hence, M is a supersolution and thus u < M a.e. on ^ by Proposition 16.71 which gives Items ([T]) 
and ([2]). 

Items ([3]) and (jH) follow from Items ([T]) and ([2]) by noting that —u is a subsolution for —F and 
— g if li is a supersolution for F and Item ([5]) follows by combining Items ([1]) and ([3]), while 
Item ([6]) follows by combining Items ([2]) and (j!]). □ 

The a priori estimate in Item ([6]) may be compared with [63[ Theorem 1.5.1 & 1.5.5] and [69^ 
Lemma 2.8]. 

Definition 6.11 (Weak maximum principle property for a linear operator). Let A : H'^{&, tt)) — )■ 
L^(^, id) be a linear operator and let S ^ d&. We say that yl obeys the weak maximum principle 
property on ^ U S if whenever u G H'^{&, Id) 

Aii < a.e. on ^, 

li < on 3^ \ S in the sense of if^(^, tt)), 
and, if & is unbounded, ess sup^ u < oo, then 

u < a.e. on 0". 

Although the application we have in mind is where A is a linear, second-order partial differential 
operator whose principle symbol is zero on S, that situation is not assumed. The comments in 
Remarks 16.41 and 16.51 also apply to Definition 16.111 
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Definition 6.12 (Integration by parts). Let A : H^{ff,m) — )■ L^{ff,'m) be a linear operator and 
let S ^ 9^. We say tliat a bilinear map a : H^{^,Xo) x H^{ff,tt)) — ;> M is associated to A on 
1^ U S through integration by parts if 

(An,?;)i2(^^n,) = a(n,?;), \fu e H\ff,iv),v e H}f{^ UJ:,m). (6.10) 

Examples of sufficient conditions for (j6.10p to hold for bilinear maps, o, and operators, A, are 
described in ^8.31 The following lemma, whose proof is clear, relates the weak maximum principle 
property on U S in Definition 16.31 to that in Definition 16.111 

Lemma 6.13 (Relationship between weak maximum principle properties). Let A : H'^{ff,w) — )■ 
L^(i^, tt)) be a linear operator, let S ^ dff, and let a : H^{ff,tt3) x i/^(^, ro) R be a bilinear 
map which is associated to A on L) T, through integration by parts. Then A obeys the weak 
maximum principle property on ^ U S if and only if the same is true for a. 

We have the following analogue of Definition 16.61 

Definition 6.14 (Strong solution, subsolution, and supersolution). Suppose / G L^(^, ro) and 
g E if^(^, ro). We define u £ H'^{&, Vo) to be a strong subsolution if 

Au< f a.e. on (6.11) 
u < g on 9^ \ S in the sense oiH^{&,Xo), (6.12) 

and ess sup^it < oo if ^ is unbounded. We call u G ff^(^, tt)) a strong supersolution if —u is a 
strong subsolution and a strong solution if it is both a strong subsolution and supersolution. □ 

We then have 

Proposition 6.15 (Weak maximum principle and a priori estimates for H'^{ff, ro) functions). Let 

f G L'^{ff,m), and g £ H^{^,n), and A : H'^{ff,w) L'^{ff,m) be a linear operator on 0' UJ^ 
associated to a bilinear map a on H^{i?',n}) through integration by parts. Assume A obeys the 
weak maximum principle property on ^US. Then the conclusions of Proposition [6.10\ hold for 
functions u G H'^{&,'m), provided the properties ()6.8p or (16. 9p for a are replaced by the properties 
for A that 

Al>0 a.e. on ff, (6.13) 
or that there is a constant cq > such that 

Al > Co a.e. on ff, (6.14) 
and the role of Definition \6.6\ is replaced by that of Definition \6.14\ 

Proof. When u is a strong subsolution (supersolution, solution), then it is necessarily a variational 
subsolution (supersolution, solution) using (|6.10p and so the result follows immediately from 
Proposition 16.101 □ 

Finally, we consider an application of the weak maximum principle property to unbounded 
functions on unbounded domains, by analogy with §2.31 Our Definition 16.31 of the weak maxi- 
mum principle property requires that the subsolution, u, on an unbounded domain be essentially 
bounded above. However, if a bilinear map has the weak maximum principle property for subso- 
lutions which are essentially bounded above, we may obtain an extension for subsolutions which 
instead obey a growth condition. 
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When u e C^(^) and v G C^{ff U S), we have a{u,v) = {Au,v)L^ff^ro)- If V G C^i^) and 
[A, (/?]u = —B{ipu), as in (j2.15p . and 

a{Lpu,v) = (^(/?u,t;)i2(^^„) 

and so 

a(v?'u,t;) + {Bipu,v)L2i^^^^^) = a{u,ipv)L2(^^^^^y (6.15) 

Let B : tt)) — t- L^(^, ro) be the first-order differential operator defined by (|6.15p for all 

u G H^{0',tv) and ?; G Hq{0' VJ Il,tt)). We then have the following version of Theorem I2.2UI for 
H'^{ff,ro) functions. 

Theorem 6.16 (Weak maximum principle for unbounded H^{i?,Xo) functions on unbounded 
domains). Let ^ Mf^ be an unbounded domain and let a be a bilinear map on H^{^,Xo). 
Assume that the bilinear map, 

a{u,v) := a{u,v) + {Bu,v)l2(^^^„), \fu,v e H^{ff,tv), (6.16) 

obeys the weak maximum principle property on U S, for some S ^ dff, for functions u G 
H^{ff,to) which are bounded above. Then, a has the weak maximum principle property on ffWE 
for functions u G H^{i^,tv) obeying the growth condition (j2.17p a.e. on G. 

Proof. Since o(ii, v) < 0, the identities (|6.15p and (j6.16p imply that a{ipu, v) < for all v G 
iJo(^ U Sjtt)), f > a.e. on ^. The conclusion now follows, just as in the proof of Theorem 
[Al □ 

Corollary 6.17 (Weak maximum principle for unbounded H'^{&,Vo) functions on unbounded 
domains). Let be an unbounded domain and let A : iif^(^, ro) — ?• L^(^, tr) be a linear 

operator on ^ U S, for some S ^ 80" , associated to a bilinear map a on H^[0',xo) through 
integration by parts. Assume that A : iJ^(^, tr) — )• L^(^, ro) in (j2.16p obeys the weak maximum 
principle property on ^ U S for functions u G H'^{0',Vo) which are bounded above. Then, A has 
the weak maximum principle property on U for functions u G i7^(^, ro) obeying the growth 
condition (j2.17p a.e. on & . 

Proof. When u G H'^{&, Vo) is a strong subsolution, then it is necessarily a variational subsolution 
since a{u,v) = (^w, ?;)^2(^ for all v G LL^^i)' U S,tt)) by (j8.3ip and the result follows from 
Theorem 16.161 □ 

7. Applications of the weak maximum principle property to variational 

inequalities 

We prove uniqueness for solutions to variational inequalities defined by bilinear maps a on 
iJ^(^, ro) and obstacle problems defined by bounded, linear operators A from H'^{ff,w) to 
L'^{ff, tt)), when a or A obey the weak maximum principle property (Definition 16. 31 or 16. IH respec- 
tively) . Applications to variational inequalities are much simpler than in the case of the obstacle 
problems considered in ^ and include, in §7.11 a comparison principle for supersolutions and 
uniqueness for solutions to variational inequalities (Theorem 17. 2p and, in §7.21 the corresponding 
results for i/^(^, tt)) solutions to the obstacle problem (Theorem 17. 7p . In §7.3^ we develop a 
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priori estimates (Proposition 17.9] ) for H^{ff,to) supersolutions and solutions to variational in- 
equalities and then the corresponding results for (Proposition I7.10|) for It)) supersolutions 
and solutions to the obstacle problem. 

7.1. Comparison principle for supersolutions and uniqueness for solutions to 
variational inequalities. We first recall our analogue [18] of the standard definition [9l |39l HSl 
[65l [69] of a solution and supersolution to a variational inequality. 

Definition 7.1 (Solution and supersolution to a variational inequality). Given a source functional 
F E H~^{ff U S,tt)), boundary data function g £ H^{0',n}), and an obstacle function ip S 
H'^{ff, to) such that ■0 < 5 on \ S in the sense of H^{0', to), that is, 

{i;-g)+ eH^i&uJ:,tv), (7.1) 

we say that a function u G if^(^, tt)) is a solution to a variational inequality with Dirichlet 
boundary condition along 80" \ S if u = <^ on 5^ \ S in the sense of H^{ff, ro), that is, 

u- g e i?o(i^U S,tt)), 

and 

u > ip a.e. on 0' and a(u, v — u) > F(v — u), 

(7.2) 

e i7^(^, tt)) with u > V a.e. on 1^ and - g G i?o (^ U S, m), 

and, if & is unbounded, ess sup^ \u\ < oo. We call a function u G H^{i^,tv) a supersolution if 
u > g on \ S in the sense of H^{0', ro), that is, 

(u-5)- Gi^o(^UE,tt)), 

and 

u> ih a.e. on ^ and a(n, w) > F(w), 

(7 3) 

G iJo(<^ U S, ro) with u; > a.e. on ^, 
and, if & is unbounded, ess inf^u > — oo. 

If is a solution in the sense of Definition 17.11 then we see that it is also a supersolution by 
writing v = u+w and observing that v > u > a.e. on and v = g on d0'\Ti if u; G i^g (^US, to) 
and li) > a.e. on 0". 

We now assume that the bilinear map a has the explicit form, 

a{u,v) := {a'^u^^v^^ +dHvx^ -b'u^^v + cuv)tt)dx, Vu, t; G C^(^), (7.4) 

for some to as in Definition 16.11 (not necessarily coinciding with the roi in the Definition 16.21 of 
H^{0, to)) and where the coefficients a*-' , , 6*, and c are Borel measurable functions on ^ C M'^. 
We require, in addition, that the coefficient matrix, a = (a*-'), obeys a version of ()2.ip for Borel 
measurable functions, namely 

{av,v) > a.e. on 0, Vr/ G R'^. (7.5) 

We then have the 

Theorem 7.2 (Comparison principle for supersolutions and uniqueness for solutions to the vari- 
ational inequality). Let F G H~^{ff\J and g G H^{0,tv), and ip G H^{0,tti), and a be 
a bilinear map of the form (I7.4p on H^{ff,to) obeying the weak maximum principle property on 
^ U S, for some T, ^ Off, and, in addition, that 

di = onff, l<j<d. (7.6) 
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Suppose ui,U2 G ((^ U S, It)) are a solution and supersolution, respectively, to the associated 
variational inequality (Definition \7.1\ ). Then U2 > ui a.e. on G and if U2 is also a solution, then 
U2 = ui a.e. on G . 

We first consider tlie simpler case where the solution, supersolution, and obstacle function are 
continuous on G ^ in which case we can adapt the proofs of |391 Theorems 1.3.3 & 1.3.4]. 



Proof of Theorem \ 7. S\ for continuous functions. For this special case, the hypothesis ()7.6p has no 
role and is not required. By replacing u hy u — g, and v hy v — g, and F hy F — a{g, •), and V' 
hy ip — g, if needed in Definition 17.11 we may assume without loss of generality that g = on G. 
We assume in addition to the hypotheses of Theorem 17.21 that ui,U2,'4' S C{G). Since U2 is a 
supersolution, then 

a{u2,w)> F{w), G iJo(^U S,rD), w > a.e. on ^, 

while ui is a solution and thus obeys 

a{ui,v - ui) > F{v - ui), Vf G iJo(^U S,rD), v > ip a.fi. on G . 

Choose V = u\ — w, where w G Hq{G U S, ro) and < w < ui — iJj a-e. on G, so that v > ip 
a.e. on G. li ui = a.e. on G, then we can reverse the roles of ui,U2 and Ui = a.e. on 
G, i = 1, 2, we would be done, so we can assume without loss of generality that the open subset 
Gi := {P £ G : Ui{P) > V'(-P)} C ^ is non-empty for i = 1 or 2. We assume that Gi / and 
note that 

a{ui,-w)>F{-w), G iJo(^U S,tt)), Q < w < m - a.e. on u. 

Thus, 

o(u2 - > 0, G -ffo(^U i;,rD), < u; < ui - a.e. on ^, 

and in particular, noting that Gi = {ui > ip} is non-empty and open and ui — ip £ C{Gi), 

a{u2 — ui,w) G -ffo(Gi U S, ro), Q < w < ui — ip a.e. on Gi. 

By the weak maximum principle property, we obtain U2 > ui on Gi <Z G and thus U2 > ui on 
^, since U2 > "p = ui on G \ Gi. 

If ui , U2 are two solutions to (j7.2p , then the fact that U2 is necessarily also a supersolution yields 
U2 > ui a.e. on G and reversing the roles of ui,U2 yields ui > U2 a.e. on G. This completes the 
proof of uniqueness for continuous solutions. □ 

For the general case, we shall adapt the proof of |69l Theorem 4.27]. 



Proof of Theorem \ 7. ^ As in the preceding case, we may assume without loss of generality that 
g = 0. Suppose that ui is a solution and U2 is a supersolution and set u := {ui — U2)~^, so 
ui Au2 = ui — {ui — U2)~^ = ui — u. Observe that u G H^{G, ro) and ui — U2 < on dG \ S and 
thus n = (ui — ^2)"'" = on dG \ S, in the sense of H^{G, ro), and hence u G Hq{G, tr). 

We claim that a(u, vu) < for all vu G Hq{G U S, ro) with it; > a.e. on G. If not, there exists 
w G C(J(^ U S), < w < 1, such that 

a(M,u;)>0. (7.7) 

For e > 0, define 

, uw 
w — 



u + e 



Observe that G U S, ro) and 

i 

'^^ ~ u + £ ' {u + e) 



uwx eux w 

< = ^ + 7^r^, l<j<d. (7.5 
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Choose = ui — ew'^ and observe that G //q (i^ U S, ro) and 

> til — ■UTU > ui — u = ui A ti2 > a.e. on ff. 
Because ui is a solution to the variational inequality (17.2p . we have 

a(ni,-y - ui) > F{v - m), Vt; G i7o((^U !!,«)) with v > ip a.e. on i^, 
and thus, choosing v = and using — ui = —ew^, we obtain 

o(-ni,u;^) > -F{w''). 
But ti2 is a supersolution to the variational inequality (17. 2p . thus 

a(n2,w^) > Fiw""). 
Adding the preceding two inequalities yields 

a{u2 — ui,w^) > 0. 
Since u = {ui — U2)~^ and = uw/{u + e), we obtain 

a{u,w^) < 0. 
Using ()7.8p , the expression (I7.4p for a yields 

a(n, w^) = J [al'^UxiWl.. + d^iiw^. — h^iixiW^ + cuw'^^ tt) dx 

= / — - — (a^^Ux Wx + d^uwx — Uux w + cuw) ro dx 

f XV 

(u + e)2 ^ ^ ^ 

=: +e/2(£)- 

Consequently, by combining the preceding inequality and identity, we obtain 

+e/2(£) < 0, Ve>0. (7.9) 

By hypothesis (|7.6p . we have {d^ ) = on and so the non-negative characteristic form condition 
(17. 5p implies that 

f w 

hi^) = / 7^— — r^a'^Ux,Ux,todx > 0. 
Jff [u + ey 

Therefore, Ii{e) < —el2{e) < for all e > and because 

Y\mli{e) = a{u,w), (7.10) 

we obtain a{u, w) < 0, contradicting our assumption that a(n, w) > 0. 

Thus, a('u,w) < for all w S HQ^ff U S,lt)) with tn > a.e. on 0" and so the weak maximum 
principle property implies that u < a.e. on i^, that is (ui — ^2)"^ = a.e. on and so 
ui < U2 a.e. on This completes the proof that if ui is a solution and U2 a supersolution to 
the variational inequality (|7.2p . then U2 > ui a.e. on i^. 

If ui,U2 are two solutions to (I7.2p . then ui = ^2 a.e. on just as before. This completes the 
proof. □ 

If we strengthen the condition (17. Sp . we can allow non-zero coefficients (d^) in (17. 4p for the 
statement and proof of Theorem 17.21 We constrain the behavior of the coefficients a, d in the 
expression (|7.4p for the bilinear map 0, near finite portions of dff as well as spatial infinity, with 
the aid of the 
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Definition 7.3 (Degeneracy coefficient). We call "d a degeneracy coefficient if 

1} e Cioc(^) and > on ff. (7.11) 

For our generalization of Theorem 17.2^ we now require that the coefficients (a*-'), (d^) in ()7.4p 
obey 

(ar/,r/) > i^lryp a.e. on 1^, (7.12) 
(ar?,r?) < K'd\r]\^ a.e. on ^, (7.13) 
\{d,7])\ < K^\r]\ a.e. on ^, (7.14) 

for all G M^. When the coefficients, (a*-'), in ()7.4p are only Borel measurable on ^, our previous 
characterization (j2.2p of the degeneracy locus is not meaningful and instead we have the 



Definition 7.4 (Characterization of the degeneracy locus for a bilinear map). For a bilinear form 
as in (|7.4p . with measurable coefficients (a*-'), the degeneracy locus, S ^ for is given by 

S = int ({x G : i9(x) = 0}) . (7.15) 

Naturally, the definitions 1^2^ and (17^5]) coincide when (a*^) G Cioc(^; M'^^'^). 

Theorem 7.5 (Comparison principle for supersolutions and uniqueness for solutions to the vari- 
ational inequality). Assume the hypotheses of Theorem \ 7.S\ but now require that either the co- 
efficients (d^) in the expression (17. 4p for the bilinear map a obey (|7.6p or that the coefficients 
(a*J),(dJ) in 1^ obey (i7J2]l . fTTSD . and fM) and that 

i9 G Li(^,tt)). (7.16) 



Then the conclusions of Theorem 7.2 continue to hold. 

Proof. We need only show that the assumption (|7.7p still leads to a contradiction using the 
inequality (j7.9p even when (d^) is non-zero. By (j7.7p . (|7.9p . and (j7.10p . we obtain 

limsupe/2(e) < — lim/i(e) = —a{u,w) < 0, (7.17) 
and so there is an eo > such that ehi^) < and hence hi^) < 0, Ve G (0,eo]- Define 

Gs{u) := — : a.e. on ff. 

Then, the fact that hie) < 0, Ve G (0, Eq], yields 

f f II )\ T~) ? / 1 ^ f TV u u 

/ ^Ge(u?K>dx= / 1?-^ ^tt)(ix< / a^^ -r^^tt) da; fbv fTUh ) 

Jff Jff [u + e)2 - Jff {u + e)2 

wuu. 

d^ - 

/ ( 



< - / c^'^T^^ TTTttifix (since 12(e) < 0) 



<K \ ^ 



wdx = K [ i)Geiu)w dx (by (mill ) 



U + £ 

1/2 / /■ \ 1/2 



<K(^J^tt>dx^ (^J ^Geiuftvdx 



1/2 

Co ( I^Geiuftvdx 
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1 /2 

where Co K [f^, ^tv dx) . But then 

||V^G,(n)||i2(^,„) <Co, VeG(0,eo], 

and so, Ve G (0, Eq], 

f w ■ ■ 

<K f ??-r^^ (\Du\^ + \Du\u) lit) dx (by (I713D and (ITliD ) 

(n + e)^ ' 

<K / -d--^ -7rXvdx + K\ I i}—^ -^rtxidx] / i}— -Trtvdx 



1/2 



<{K + l)Cl 

Therefore, |e/2(e)| < {K + '^)Cq£ and so hnie^o £-^2(e) = 0, contradicting (j7.17p . □ 

7.2. Comparison principle for supersolutions and uniqueness for solutions to 
obstacle problems. We recall our analogues of the standard definition [9l [39} HGj [65| [69] of a 
strong solution to an obstacle problem defined in |18j . 



Definition 7.6 (Strong solution and supersolution to an obstacle problem). Given functions 
/ G L2(^, ro), g E i?i(^, ro), and -0 e if^(^,rD) obeying (fTTT]) . we call -u G ii"2(^, tr) a strong 
solution to an obstacle problem for a linear operator, ^ : H'^{0',tv) — > ^^(1^,03), with Dirichlet 
boundary condition along dff \ S, for some S ^ dff, if 

min{Au - /, n - ^} = a.e. on ^, (7.18) 
ti-g G F(J(^US,ro), (7.19) 

that is, u = g on 9^ \ S in the sense of H^{ff, n), and, if is unbounded, ess sup^ \u\ < 00. We 
call u G -ff^((^, ro) a strong supersolution if 

min{An - /, n - ^} > a.e. on 1^, (7.20) 
(n-5)- GFi(^US,ro), (7.21) 

that is, n > 5 on dff \ S in the sense of H^{i?, tu), and, if i?' is unbounded, ess inf^ u > —00. 

Theorem 7.7 (Comparison principle and uniqueness for solutions to the obstacle problem). 
Let f G L^{ff,w), and g G H'^{^,m), and G H^{^,Xo), and A : H^{ff,n) L'^{ff,Xo) he a 
linear operator on ffUT, associated to a bilinear map a of the form ()7.4p on H^{0',)xi) through 
integration by parts, for some S ^ dff. Assume A obeys the weak maximum principle property 
on ffUT, (Definition \ 6.11\) and that the coefficients (a*-'), (d-^) in the expression ()7.4p obey the 
hypotheses of Theorem \ 7. 5\ Suppose ui G H'^{ff,m) is a strong solution and U2 G H'^{ff,Xv) is a 
strong supersolution to the obstacle problem in the sense of Definition \ 7. ff| Then U2 > ui a.e. on 
G and ifui,U2 are solutions, then U2 = ui a.e on G . 

When the solution, u\,U2 G H^{G, to), are also assumed to be continuous on G, one can adapt 
the proof of \39\ Theorem 1.3.4] to prove Theorem 17.71 
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Proof of Theorem \ 7. 7| for continuous solutions. We assume in addition to the hypotheses of The- 
orem O that ui,U2 G C(^). Let = {P G ^ : U2{P) > ui{P)} C ff. (If ^ is empty, then 
ui < U2 on 0" and if ui,U2 are both solutions, we may reverse the roles of ui and U2 to give 
ui = U2 on Since M2 > ui > ■;/' on then Au2 = / a.e. on and Aui > f a.e. on 
implies that A{u2 — ui) < a.e. on . Since ui — U2 = on n 9^, the weak maximum 
principle property implies that U2 < ui on a contradiction. □ 

For the general case, we shall use the following analogue of the equivalence in [9j Equation 
(3.1.20)], whose proof is identical to the proof of the corresponding equivalence lemma in |18] . 

Lemma 7.8 (Equivalence of variational and strong (super-)solutions). Assume the hypotheses 
for the operator, A, in Theorem \7. 7\ Let f , and g, and ip be as in Definition \7.6\ and suppose 
u G H'^{ff, Vo). Then u is a (super-) solution to the variational inequality in Definition \7.1\ if and 
only if u is a (super-) solution to the obstacle problem in Definition \ 7. 6) . 

Proof of Theorem \ 7. 7\ Since ui,U2 are a strong solution and supersolution in the sense of Def- 
inition 17. 6^ then they are necessarily a solution and supersolution to the variational inequality 
in the sense of Definition 17.11 by Lemma 17. 8|, and so ui < U2 a.e. on by Theorem 17. 5| while 
ui = U2 a.e. on if both ui and U2 are strong solutions. □ 

7.3. A priori estimates for solutions and supersolutions to variational inequalities. We 

now state maximum principles and corresponding a priori estimates which extend those of non- 
coercive variational inequalities defined by uniformly elliptic partial differential operators, such 
as [Ml Theorems 4.5.4 & 4.7.4] (see also [Ml Theorem 4.5.1 k Corollary 4.5.2]). If n G H^{^, ro) 
is a supersolution to the variational inequality in Definition 17. H then u necessarily obeys the 
inequality, 

u < ip a.e. on ff, 

but the weak maximum principle yields additional a priori estimates for u, as described below. 

Proposition 7.9 (A priori estimates for supersolutions and solutions to variational inequalities). 
Let F = {f,f\..., f^) G H-^{& U S, ro), and g G H^{&, Xo), and ^ G H^{ff, Xo), and a be a 
bilinear map of the form (|7.4p on H^{iff,xo) obeying the hypotheses of Theorem \ 7. 5[ for some 
Tj ^ dff, and assume that (16. Sp holds. Suppose that u G H^{ff,to) is a supersolution to the 
associated variational inequality (Definition \7.1^ . 

(1) IfF>0, then 

M > A ess inf g a. e. on G . 

(2) If F = (/, 0, . . . , 0) and f has arbitrary sign and a obeys (j6.9p . then 

u > A — ess inf / A ess inf g on G . 

Co d0\T, 

(3) If F < and u is a solution for F and g and i{j (Definition \7.1^ , then 

< V ess sup^i V ess sup^/^ a.e. on G . 

(4) If F = (/, 0, . . . , 0) and f has arbitrary sign, u is a solution for F and g and if), and a 
obeys (|5.9p . then 

ti < V — ess sup / V ess sup g V ess sup il) on G . 

Co e dff\T: ff 
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(5) If ui and U2 are solutions, respectively, for Fi > F2 and tpi > ^2 o-e. on G , and g\ > §2 
on dG \'S in the sense of H^{G, m), then 

ui > U2 a.e. on G . 

(6) If u\ and U2 are solutions, respectively, for Fi > F2 and ipi > ^2 o-e. on G, and gi > g2 
on dG \ S m the sense of H^{G, tv), and a obeys (I6.9p . and Fk = {fk,0, • • • , 0), k = 1,2, 
then 

- U2\\l°°{6') < — ll/i - /2||l°o(^) V \\gi - 52||l°°(9^\S) V 11^1 - V'2||l°o(^), 
Co 

while if fi = f2 a.e. on G and a obeys ()6.8p . then 

\\ui - U2\\L°^(ff) < \\gi - 52||l°o{9<?\S) V Wlpl - V'2||l°°(<5')- 

The terms ess supg^x^-^g, and ess infg^^^^; o.nd ||5'||l°°(9^\s)) o.^d \\gi — 52||loo(9^\s) iiT- 
preceding items are omitted when S = dG. 

Proof. Consider Items ([1]) and ([2|). Since n is a supersolution to the variational inequality in 
Definition 17.11 then it is also a supersolution to the variational equation in Definition 16.61 (where 
ip plays no role) and so Items ([T]) and ([2]) here just restate Items ^ and in Proposition 16. 101 
Consider Items ^ and (g]). When / < a.e. on G, let 

M := V ess sup g V ess sup ip, 

while if / has arbitrary sign and a obeys (|6.9p . let 

M := V — ess sup / V ess sup g V ess sup 

Co e 

We may assume without loss of generality that M < 00. Then M > ip a.e. on G and M > g on 
dG \ E in the sense of H^{G,n), while for ah w G H^{G\Jj:,w) with w>0 a.e. on G, we have 

a{M,w) = (cM,u;)i2(^^^„) > > (/, w)j^2(^^„), 

when / < a.e. on G and 

a{M,w) > (coM,'u;)i2(^^„) > (/, w;)l2(^^„), 

when / has arbitrary sign. Hence, M is a supersolution and so Theorem 17.51 implies that u < M 
a.e. on G, which establishes Items ([3D and dH). 

The proofs of Items 1^ and ([U]) here are very similar to the proofs of the corresponding Items 
(0) and ([HD in Proposition 13. 5| except that appeals to Proposition 13.31 are replaced by appeals to 
Theorem 17.51 □ 

Note that an L°° comparison estimate for solutions ui,U2 corresponding to /, g, and obstacles 
V'i,V'2 is provided by [65l Theorem 4.7.4] and [H Theorem 3.1.10]. 

Proposition 7.10 (Weak maximum principle and a priori estimates for strong solutions to 
obstacle problems). Let f G L'^{G,{o), and g G II^{G,{o), and A : H'^{G,w) L'^{G,K)) be a 
linear operator on GUTi associated to a bilinear map a on H^{G, ro) through integration by parts. 
Assume A obeys the weak maximum principle property on GuT,. Suppose that u G H'^{G,tv). 
Then the conclusions of Proposition \ 7.9\ hold, provided the properties ()6.8p or (16. 9p for a are 
replaced by the properties (|6.13p or (j6.14p for A, the role of Definition \7.1\ is replaced by that of 
Definition \ 7. and supersolution is replaced by solution. 
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Proof. When u is a strong (super-)solution, then it is necessarily a variational (super-) solution 
using (|6.10p and so the result follows immediately from Proposition I7.9[ □ 

8. Weak maximum principle for functions in Sobolev spaces 

Having considered applications of the weak maximum principle property (Definition 16. 3p to 
variational equations in ^and variational inequalities in ^ we now establish conditions under 
which the bilinear map a on H^(ff,K)) or the associated differential operator A : H'^{0',K)) — >■ 
L^(^, tt)) has the weak maximum principle property on ^ U S. We begin in ^8.11 with some 
technical preliminaries. In ^8.21 we prove a weak maximum principle for o when i7^(^, tr) is 
defined by power weights and ^ is a bounded subdomain of the upper half-space, EI C M"^ 
(Theorem 18. 8p : modulo a suitable weighted Sobolev inequality (Hvpothesis 18. 9p . we then prove 
a more widely applicable weak maximum principle for a when ff^(^, tt)) is defined by general 
weights and ^ is a bounded subdomain of (Theorem l8.1ip . In ^8.31 we describe an integration 
by parts formula relating, under suitable conditions on the weights and coefficients, the bilinear 
map, 0, on i?^(^, tt)) and the associated operator. A, on if^(^, ro). Finally, in ^8.41 we prove 
a weak maximum principle for bounded H^{0',Xo) functions on unbounded domains (Theorem 

EM- 

Our weak maximum principle differs in several aspects from |63t Theorems 1.5.1 &: 1.5.5], which 
again may appear subtle at first glance but which are still important for applications: 

(1) We use weighted Sobolev spaces adapted to the coefficients of the first and second-order 
derivatives in A and so those conditions are weaker than those of |63l Theorems 1.5.1 &: 
1.5.5] and permit applications to operators such as those of Examples 11.21 and 11.41 which 
are not covered by |631 Theorems 1.5.1 & 1.5.5]; 

(2) The subdomain ^ C is allowed to be unbounded. 

The differences between results obtainable from our maximum and comparison principles and 
those of Fichera, Oleinik, and Radkevic are described in more detail in gO using the example of 
the Heston operator (Example II. 2p . 

Remark 8.1 (Strong maximum principle for bilinear maps on i7^((^. It))). We shall not consider 
in this article an analogue of the strong maximum principle for a on i7^(^, ro) corresponding 
to |40t Theorem 8.19], as such a result requires a corresponding version of the weak Harnack 
inequality Theorem 8.18]. In the case of the Heston operator (Example ll.2p and its associated 
bilinear map, both a weak and a strong Harnack inequality (analogous to Theorem 8.20]) are 
proved in [31] and we expect that those Harnack inequalities extend to the setting considered in 
this section. 

8.1. Preliminaries. We let Xo be as in Definition 16.11 and i? be as in Definition 17.31 and obey 
()7.16p and choose 

M\m(e,ro)-= I {^\Du\^ + {l + nu?)^dx, (8.1) 

so that roi^o = (1 + ^)tt' aiid tui^i = "dtv in the Definition 16.21 of H^{0', ro), while tt)o,o = tt) in the 
definition of L^(^, Id). 

We shall assume that the coefficients, (o*-') and (d-'), of o in ()7.4p obey (|7.13p and (17.14p and, 
in addition, that the coefficients (6*) obey 

\{b,7])\ < K{}\r]\ a.e. on ^, V?7 G R'^, (8.2) 
\c\<K{l + 'd) a.eon^, (8.3) 
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for some positive constant, K. It is easy to check that a obeys the continuity estimate, 

a{u,v) < Ci\\u\\Hi{ff,n)M\m{ff,n), yu,v e H'^{ff,Ki), (8.4) 

for some positive constant, Ci, when the coefficients, (a*-'), (d^), and (6*), of o in (j7.4p obey (|7.13p . 
(HUD, dO]), and (HSI), in which case Ci = Ci{K) in (18^ . 

For 1 < p < oo, we let L^'(^, id) denote the Banach space of Borel-measurable functions, u, on 

such that 

Mlp{6\w)-= \u\Pxodx<oo. (8.5) 

J & 

The bilinear map a obeys a Gdrding inequality^ 

a{u,u)>C2\\u\\l,^0^^^-Cs\\{l + ^Y/MLHfi», VnGFi(^,tr), (8.6) 

if, in addition, we require that the coefficient matrix, a = (a*-^), obeys (j7.12p . in which case 
C2 = C2iK) and C3 = C3{K). 

We now specializj^ the Hilbert space i/^ (€?,«)) in Definition 16.21 to be the completion of the 
vector space C^{d) with respect to the norm (jS.ip . Given S ^ 3^, we let H^{ffU S, ro) be as 
in Definition 16.21 

8.2. functions on bounded and unbounded domains. We first consider a special case 
of our desired maximum principle, analogous to [40 ^ Theorem 8.1]. 

Theorem 8.2 (Weak maximum principle for H^{0',tv) functions). Let he a bounded 

domain and let S ^ di? . Assume that the coefficients of in (17. 4p are Borel-measurable on 6 , 
obey ^CfXh . (l7J3]l . (17141) . and require that tt) and-d obey 

infi^ro > and supi?tt) < 00. (8.7) 

Ifue H^{^,Xo) obeys ([621) with / = 0, then 

ess sup n < V ess sup u. 

Moreover, a has the weak maximum principle property on 1^ U S in i/ie sense of Definition \6.3X 

Remark 8.3 (Boundedness requirement on ff). In applications, the functions ■!? and to would 
not normally obey the upper bound in (j8.7p unless were bounded and the lower bound in 
()8.7p unless a = (a*-') were uniformly elliptic on ff. Moreover, boundedness of is implicitly 
used in the proof of Theorem 18.21 in inequalities involving the Lebesgue measures of the supports 
of functions and their gradients. However, unlike the proof of \A0\ Theorem 8.1], our proof of 
Theorem 18.21 avoids the use of the Poincare inequality, thanks to a nice observation of Camelia 
Pop, and hence a requirement that ^ is bounded originating from the usual statements of the 
Poincare inequality (for example, j28|. Theorem 5.6.3]). 

Proof. We essentially follow the proof of |40| Theorem 8.1], but include the details here for 
later reference in our proof of Theorem ESI If u G H^{ff,m) and v £ H^{ff U then 
uv € Wq'^^i^ U S, tv) and D{uv) = vDu + uDv by analogy with \AQ\ Problem 7.4], recalling that 



We prove an analogue of the Meyers- Serrin theorem [3] for unweighted Sobolev spaces in the context of certain 
weighted Sobolev spaces in [18| . 
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H\ff,tt>) = VFo'^(^,tD) and H^{ffU^,Xo) = VFq '^(^ U S, w). By the definition ^ of the 
bilinear map, a{u,v), we obtain 

/ (a^-'uxiVx. — {d^ + V')ux^v) Xvdx < — / (d^{uv)x, + cuv) ro dx < 0, 
Jff Je 

for all V G i/Q(€?'US, xo) such that f > and itf > a.e. on ^, where to obtain the last inequality 
we use (|6.8p . Therefore, 

/ a^Jna^^Vx.rodx < / ((i* + 6*)na;,utrdx < K / ^\Du\v^dx fbv (f7l4]l and ^JS^ \. 
Jg Je Je 



Denote 



/ := V ess sup it > 0, 



and recall our convention (I6.3P that V ess supg^^j] u = when dG \ Yj = 0. We may assume 
without loss of generality that / < oo, as otherwise there is nothing to prove. Suppose there exists 
a constant k such that 

/ < A; < ess supn < +oo. (8.8) 

(If no such k exists, then we are done.) Set 

^;:=(n-/c) + , (8.9) 

and observe that, because u < k on dff \ S in the sense of H^{ ff, tv) (when dff \ S is non-empty) 
and u,k e H^{ff,ro), then v e H^{ffUT.,m) with 



Dv 



Du for u > k, 
for u < k. 



Consequently, if ^ denotes the interior of suppD« C suppu, 

/ '&\Dvfxodx<K i!}\Dv\vn dx, 

since {Du,Dv) = \Dv\'^ a.e. on iff and ()7.12p gives 

/ a^^ VxiVx to dx > / '&\Dv\'^wdx. 
Je ' ^ 

Therefore, by the Cauchy-Schwartz inequality, 

\W^'^Dv\\l2^^,^„) < Jf||^V2^|U2(.^,„). (8.10) 

The inequality ()8.10p yields 

\\Dv\\L2(a^'^ < Ci||t;||£,2(^^), (8.11) 

for a positive constant 



Ci := K 



f sup^. -dXP \ 
V va.iff'dxo ) 



1/2 



which is finite by (j8.7p . Combining the inequality ()8.1ip with the Sobolev embedding Wq ' (fi/ ) — )■ 
L9('^), 2<g<ooif(i = 2and2<g< 2d/(d - 2) if d > 2 [3l Theorem 5.4 (Parts I (A & B) & 
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III)] and the fact that v G Wq''^{^) imphes, for positive constants C2, C3 depending on Ci, q, , 
that 

lb ^ C2 (||-Du||j^2(o;/) + ||w||l2(<)^)) (Sobolev embedding with q > 2]) 

< <^3lbllL2(^^) (by (EH])) 

< C3|^|^/^^^/''|b||L,(.^) (by Uni Equation (7.8)]), 
and thus, recalhng that is the interior of suppDu, 

|suppZ)i;|^/2-^/« > C3 > 0. 

The inequality is independent of k < ess sup^ti and so continues to hold when we take the 
limit k — )• ess sup^ti. Hence, we see that u attains its maximum ess sup^ti < +00 on a set, 
, of positive measure. If ess sup^u = +00 on we obtain a contradiction to the fact that 
u G L^('^), since u G L^(^, tr); if ess sup^ u < +00 then, because u is constant on , one must 
also have Du = on , contradicting the fact that 

1^ n suppDn| = I suppL't'l > 0. 

This completes the proof. □ 

We now suppose 0' cM and recall the following 

Theorem 8.4 (Weighted Sobolev inequality for power weights). [471 Theorem 4.2.2] Suppose 

1 ^ P ^ q < 00, and s > —1/p, and 1 — d/p < ^ < 1 is defined by 

- = - + — ^. 8.12 

p q d 

Then there is a positive constant C = C{d,p,q,s) such that, for any u G L'?(]HI, x^) with Du G 
LP(]H,x^+";M'^), one has 

||2;>||i9(H) < C||x^+^L'u||iP(H). (8.13) 

Remark 8.5. There is a minor typographical error in the statement of |471 Theorem 4.2.2]. 

We have the following generalization of j47l Lemma 4.2.4]. 

Corollary 8.6 (Application of weighted Sobolev inequality for power weights). Suppose /3 > 
and2-d<a<2. For u e L'^{M, xJ"^) with Du G L'^{M, M'^), and q > 2 defined by 

1 1 l-a/2 , , 

and 2 < r < q, one has 
where AG [0, 1] is defined by 

1 A 1-A 

Proof. When s = (/3 - l)/2, so /? > when s > -1/2, and ^ = a/2, so 1 - (i/2 < ^ < 1 when 

2 — ci < a < 2, it follows from Theorem 18.41 that 
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Holder's inequality, in the form of [401 Equation (7.9)], gives 

II^IIl'-(h,x^1) ^ ^II^IIl2(h,x^1)II^IIl^(h,x^1)' 

when A G [0, 1] is defined by (j8.16p . Combining the preceding two inequalities yields the result. □ 

Remark 8.7 (Weighted Poincare and Sobolev inequalities). Direct weighted analogues of the 
standard Sobolev inequalities, with weights given by positive powers of the distance to a bound- 
ary portion, S ^ 9^, such as [491 Theorems 19.9 & 19.10], only hold for a very restrictive range 
of powers, even when G is bounded. Koch provides a weighted Poincare inequality on H with 
a weight similar to ours |47l Lemma 4.4.4], as well as certain weighted Sobolev inequalities [471 
Theorem 4.2.2 &: Lemma 4.2.4]. Adams [3l §6.26] provides an unweighted Poincare inequality 
which is valid on unbounded domains of "finite width" , while the Gagliardo-Nirenberg-Sobolev in- 
equality |28|. Theorem 5.6.1], a Poincare-type inequality on M"^ and the Caffarelli-Kohn-Nirenberg 
inequality, another weighted Poincare inequality on M"^ [21], are potentially useful in this context. 



With the aid of Corollary 18.61 &s noticed by Camelia Pop, we can relax the non-degeneracy 
requirement (18. 7p in the hypotheses of Theorem | 



Theorem 8.8 (Weak maximum principle for ff^(i^, tt)) functions on bounded subdomains of 
the upper half-space). Assume the hypotheses of Theorem \8.^ except that we now omit the 
requirement (18. 7p . In addition, assume G '^M is a bounded subdomain and require that there 
are constants, < c^ < 1 and < Ctx, < 1, such that 

c^x2 < 1? < c^^x;^ on ff, (8.17) 

Cmx'^d'^ <tv< CnJ^x^"^ on G, (8.18) 
where P > andO < a <2. IfueH^{G, ro) obeys ([62]) with f = 0, then 

ess sup u < V ess sup u. 

Moreover, a has the weak maximum principle property on ^ U S in the sense of Definition \6.3X 



Proof. We proceed as in the proof of Theorem 18.21 except that we now appeal to Corollary 18.61 
in place of the standard Sobolev inequality [3l Theorem 5.4]. From (|8.10p and the fact that v in 
()8.9p extends by zero outside suppv to an element of H^(M.,x^^~^) by the analogue of O Lemma 
3.22], also denoted by v, we obtain 

where Ci := {c^Cn)~^K and, as in the proof of Theorem 18.21 the set <Z & denotes the interior 
of suppDf. Then 

ll-llL^(H,.r^) < (by dm) 

^^3||^|li.(H,,,a-.)||Hlil('V,,,-i..) (bydHSD) 

^ ri {1-^)q/2|| ||A II ||1-A 

<C,smaxx\ uLo,,, fl_,Jk;" 
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for a positive constant C4 = C3 max^.^^ x^^ A)a/2 jj^^^gpgj^jg]^^ Qf ^]^g constant k in (|8.8p . where 
we apply Corollary 18.61 with r > 2 (which is possible since a < 2 and thus q> 2) and < A < 1. 
Recah from gOl Equation (7.8)] that 



I . 1^ |l/2-l/r|| II 



where we denote |^|^ := fyX^ dx for any /5 > and Borel-measurable subset ^ C H. Hence, 
the preceding inequalities give 



< r7J sunn 7/(^/2-^/'') Il7;ll^ \'9/\^^-^^^^^^~^^'''h\v\\^^^ 
<04|supp7;|^ II^IIl'-(h,x^1)I'^I/3 II^'IIl'-c^.x^I) 

-04|SUppV|^ KI/3 



and so, noting that is the interior of supp-Dw, 



Q|supp.;|f /^-^/'-)|suppI)^|(^-^)(^/^-^/'-) > 1. 
Thus, since | suppu|^ < < 00 and > 0, 

recalling that < A < 1 and r > 2. We again obtain a contradiction, after taking the limit 
k — >• ess sup^ u, and the result follows. □ 

Given a suitable weighted Sobolev inequality for functions on subdomains ^ M"' (see Corol- 
lary IB. 21 for a general class of examples for functions on W^) , it is a straightforward to generalize 
Theorem 18.81 from the case where and tt) obey (j8.17p and (j8.18p . 

Hypothesis 8.9 (Weighted Sobolev inequality). Given a subdomain ^ M'^, constants 1 < p < 
q < 00, and functions It) S C(i^) such that t} > and ro > on 1^, there is a positive constant 
C = C{p,q,'&,K)) such that, for any u € L'?(^,tt)) with Du G LP(^, i?lt); M'^), one has 

II^IIl9(^,w) ^ C\\Du\\ip(^ff^^ny (8.20) 

Corollary 8.10 (Application of weighted Sobolev inequality). Assume Hypothesis \8. 91 holds with 
p = 2. For u £ L'^{^,w) with Du E L'^{ff,'&to;R'^) and 2<r<q, one has 

\Ml'-{(^,w) < C|l^llL2(<:?,re)ll-^^lli2{^^,i»m)' (8-21) 
where A G [0, 1] is defined by (|8.16p . 

Proof. The proof is similar to that of Corollarv 18.61 Inequality (|8.20p . with 2 < q < 00, gives 

< C\\Du\\i2(^ff^^„y 
Holder's inequality, in the form of [IHl Equation (7.9)], yields 

when A € [0, 1] is defined by (I8.16p . Combining the preceding two inequalities yields the result. □ 
We have the following generalization of Theorem 18.21 
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Theorem 8.11 (Weak maximum principle for H^{if,tx}) functions on bounded domains and 
general weights). Assume the hypotheses of Theorem \8.S\. except that we allow & <^W^ to he any 
bounded subdomain with S ^ dff and allow 'd,n> to be any functions obeying Hupothesis \8. 9\ for 
p = 2 and some 2 < q < oo. If u £ i?^(^, tt)) obeys (j6.6p with / = 0, then 

ess sup n < V ess sup u. 

a0\T. 

Moreover, a has the weak maximum principle property on S m the sense of Definition 1 6. 3l 

Proof. The proof is almost identical to that of Theorems 18.21 and 18.81 except for a few minor 
changes which we indicate here. In place of ()8.19p . we note that (I8.10p may be written as 

l|-D^llL2('?/,tfro) < K\\v\\L2(a//^^^y (8.22) 

We now proceed as in the proof of Theorem 18.81 but apply (j8.2ip in place of (j8.15p and choose 
C4 = C'a maxj-g^ i^^^""^)/^ instead of Camax^-g^x^^ A)a/2^ |--| 

8.3. Integration by parts formula. Before proceeding to consider when the weak maximum 
principle holds on unbounded domains, we shall need to introduce an integration by parts formula 
and so, to accomplish this, we must impose additional conditions on the coefficients of a beyond 
those stated in §8.1[ For now, we shall require that a = (a*-') be continuous on 0' (not merely 
Borel measurable on ff) but shortly strengthen this condition further. Let A be the partial 
differential operator given in equivalent divergence and non-divergence forms by 

Au := - {a'^u^^ + d^u)^^ - [V + (log tt))^.^a*^') n^, + cu 

U + a^, +d^ + (log tt)),, a'^) n,, + (c - 4^ - (log m),.^. d^) u (8.23) 

= -a'-^Ux^x.j - b'ux^ + cu, 

with 

U := U + 4^. + (log tt))^^.a'^ l<i<d, (8.24) 
c:=c-4.. -(logw),^,d^ (8.25) 

and where we now impose the additional regularity requirements, 

«u g c^,\i^) n aoc(^), l<i,3<d, (8.26) 

d^ E CO'H^) n aoc(^), 1 < J < (8.27) 

logtD G C^'H^). (8.28) 
Provided we also require that id obey 

i9tt) G Cioc(^), (8.29) 

i9tt) = on S, (8.30) 



and that 1^ is a domain for which the divergence theorem holds, then integration by parts in (17. 4p 
yields the integration by parts relation (j6.10p when u G C^{i?) and v G C^{ff US), that is, 

a{u,v) = {Au,v)l2(^^^^„), (8.31) 
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Since 



a(u, v) = {{a'^Uxi + d?u) v^. - b'^Ux^v + cuv) W dx (by ([73|)) 

= j {~ + d^u)^^ - UUxi + CV^ VXO dx 

[a^^Ux + d-'u) vtVx dx — nj (^a^^Ux + d^u) vXv ds 

+ — d^x- — (log ^)xj d?^ dx 

= {Au,v)L2(ff^^), {hy^M) 

where n is the mwarfi-pomtmg normal vector field and the integral over di? is zero since w = 
on df \ S and tt), 1? obey ^^?M and (f830D . and the coefficients (a*J), (d^) obey (I7J2]1 and (I7J3]1 
on together with and ([8:271) . 

Remark 8.12 (Relaxing the conditions on di^). More generally, if the divergence theorem is not 
assumed to hold for i^, then (j8.3ip still holds under slightly stronger regularity assumptions on 
the coefficients, a*-^, and weight, tt), near S. 

Clearly, when the coefficients (o, 6, c) of A obey ()7.13p and 

|6| < K{1 + ^) a.e. on ^, (8.32) 
|c| < K{1 + ^) a.e. on ^, (8.33) 

there is a positive constant, C4 = C4(i^), such that 

\\M\L^{eM < C4u\\HHff,v), Vn G C^{^), (8.34) 

where we set 

\Mh(ff,vo) ■= j {Ad'^u\^ + (1 + ^'^) [\Du\^ + I^^P)) W dx, (8.35) 

so that rD2,o = ^2,1 = (1 + ^?)^rt' and tt)2,2 = ^?^rti in the Definition [O of //^((^, ro). 

We specialize the Hilbert space i/^(^, ro) in Definition 16.21 to be the completion of the vector 
space C^{ff) with respect to the norm ()8.35p and note that ()8.34p continues to hold when 
u G H'^{0', to). Furthermore, the proof of |181 Lemma 2.33] (integration by parts) adapts to show 
that ([83T]) continues to hold when u € H'^{ff,ro) and v G H^{^U'S,to). 

Remark 8.13 (Conditions on the coefficients of A and the integral weight). The bounds (|8.32p 
and (|8.33p hold, for example, if we strengthen the conditions (|8.26p . (|8.27p . and (|8.28p by requiring 
that \a%\ < K{1 + ^) a.e. on ff, and < i^(l + i9) a.e. on ff, and \(logw)xj \ < /^(l + t?) a.e. 
on ff. 

Example 8.14 (Heston operator, bilinear map, and integration by parts formula). We show in 
fCl that the coefficients of the Heston operator, A, its associated bilinear map, a, and weight 
function, tt), defined in [18], and the degeneracy coefficient, i}, obey the conditions described in 
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8.4. functions on unbounded domains. By adapting the proof of the maximum principle 
for bounded functions, Theorem 15.31 and appealing to Theorem 18. instead of Theorem 15.11 
we obtain 

Theorem 8.15 (Weak maximum principle for bounded H^{i^,ix)) functions on unbounded do- 
mains). Let G ^ he a possibly unbounded domain such that the divergence theorem holds. 
Assume the hypotheses of Theorem \8.8\ (power weights and subdomains of M) or Theorem \8.11\ 
(general weights and subdomains o/ M'^/ Assume, in addition, that c obeus^i (I2.14P a.e. on G 
and that a'^,d^,'m,'d obey IK26\i . (fOTD . <^M), (l839]l . (fS^OD . and 

/ (l + |x|2)(l + i?2)tt)(ix<oo, (8.36) 

and that ()5.7p is obeyed a.e. on G by (a, 5) in place of{a,b), where b is given by ()8.24p . Suppose 
f e L^{G, w) and sup^ f < oo. If u £ H^{G, n) obeys I^M) and ([621) with g = (when 9^ \ S 
non-empty) and 

ess supn < oo, (8.37) 

then 

ess sup n < V — ess sup /. 

ff Co & 

Moreover, a has the weak maximum principle property on GUT, in the sense of Definition 1 6. AT . 
Proof. We proceed almost exactly as in the proof of Theorem 15.31 and choose 

M := V ess sup(/ + 2Ku), 

where -fC > is the constant arising in the proof of Theorem 15.31 Our hypotheses on / and u 
imply that < M < oo. For a constant A > 0, set 

a\{ui,U2) := a(ni,M2) + A(tii, 'U2)l2(^ „,) V'Ui,n2 G i7^(i^,tt)). 

Let vq be as in dSS]) and note that vq G H'^{G,xo) by ([835]) and ([8:36]) . For (5 > 0, choose 
w G H'^{0', Xo) as in ([5.10p and observe that for all w G (^ U S, w) with i; > a.e. on G, 

a2K{w,v) = a2K{u,v) - Sa2K{vo,v) - (cq + 2K)~^a2K{M,v) (by (|5.10p) 
= o(u, v) + 2K{u, v)L2^ff^^^ - 5{{A + 2K)vo, v)L2(^ff^^^ 

- (CO + 2K)-\{A + 2K)M, t;)L2(^,„) 

(since vq,M e H'^{G,Vo) and applying ([OT]) ) 
= a{u, v) + 2K{u, v)L^()\m) - ^((^ + '^K)vo, v)L^(?,m) 

- (CO + 2K)-\ic + 2K)M, 7;)l2(^» 

(by (|2TI]1 . dO]) . and ([OD ) 

< I ess sup(/ + 2Ku)+, v) - (M, f )l2(^^„) = 0, 

where the final equality follows from the definition of M. We now apply Theorem 18.81 or 18.111 
(instead of Theorem 15. ip and the remainder of the proof is the same as that of Theorem 15.31 □ 



l^This is equivalent to (jO)) 
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Appendix A. Weak maximum principle for unbounded functions 

Theorems 12.201 and 16.161 and Coroharv 16 . 1 71 gave sufficient conditions describing when the weak 
maximum principle holds to unbounded C^(^), H^(0',tt>), or H'^{0',n) functions, respectively, 
obeying a growth condition (j2.17p defined by a function ip on the unbounded domain, ^ M.'^. 
In this section, we give examples of such growth conditions in Theorem 12.201 primarily in the 
case of the Heston operator, A. 

We first compute the coefficients of the operator B in p.lSh explicitly, when A is as in ()1.3p . 
For V G C2( 



[A, ip]v = A{{pv) — ipAv 

= {-a'^ip^^^. - b'lpx^ + cip)v+ {-a'^ {^x.Vx^ + (fx^v^^ + (fVx^Xj) - h\vx^) - (pAv 
= -a'^ {^xiVxj + ^XjVxi) - (a*Vx,x, + b'ifxi) v 
= - [a^^ + a-'*) ^XjVxi - {a'^^xiXj + &Vx,) v 
= - {a'^ + a^') ip~^(pxj {{v>'(j)xi - ^x,v) - {a'^^PxiXj + ^Vx,) v 
= - {a'^ + a^') {logip)xjMx, - {a'^^x.xj + 6Vx, - {a'^ + a^') {\ogip)x^'Px?j v, 
and therefore, since B{ipv) = —[A, ip]v by ()2.15p . we see that for all v £ C^(^), 
Bv = fvx, + fv = {a^^ + a^^) (log ip)x,Vx, 

+ (a*V" Vx,x, + b'{\og ip)x, - {a'^ + a^') (log (p)x^ (log (f)x,) v. 
Next, we give some examples of choices of functions, f. 
Example A.l (Exponential-affine growth). When if has the form 

(^(x) = e-^'^'^^, VxGM^ (A.2) 
for a fixed vector, h G M'^, and positive constant, C, the expression (jA.ip for Bv simplifies to 
Bv = - (a*^' + a^'*) hjVx, + {a^^hihj - b^hi - (a*^' + a^'*) hihj) v, 

and thus, 

Bv = -{a'^ + a^')hjVx^- {b'h, + a'^h,hj)v, v £ C^{ff). (A.3) 

Therefore, A = {A + B) is given by 

Av = -a'^Vxix^ - {b' + (a'^' + a^')hj) Vx, + (c - b'K - a'^hihj) v. (A.4) 

Hence, when ip is as in ()A.2p . it is easy to tell when A obeys the hypotheses of Theorem 12.201 
Indeed, it suffices to ensure that the coefficient, 

c := c — b^hi — a^^hihj, 

in the expression for A obeys (j2.14p for suitable h. □ 

Example A.2 (Elliptic Heston operator and exponential-affine growth). Suppose h = {hi, /12) = 
{L,N), where L > and > 0. From the identification of the coefficients, {a,b,c), for the 
elliptic Heston operator. A, in Remark 15.41 we see that 

c = c — b^hi — a^^hihj 
= r-{r-q- y/2)L - ^{6 -y)N-^ {L^ + 2paLN + a^N^) 

= I (L + 2kN -L^ - 2paLN - u'^N^) + r- k0N - (r - q)L 



(A.7) 
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Therefore, provided the coefficients obey 

L + 2nN -L^ - 2paLN - a'^N'^ > and r - kON - (r - q)L > 0, (A.5) 
we see that c obeys ()2.14p . as desired. □ 
Example A. 3 (Exponential-quadratic growth). When if has the form 

(^(x) = e-^l^l', VxGM^ (A.6) 
for some positive constant, L, the expression (jA.ip for Bv simplifies to give 

Bv = -2L {a'^ + a^') XjV^^ + {2L^a'H,j - 2Lb'xi - AL^ [a'^ + a^') XiXj) v. 
Therefore, in this case, A is given by 

Av = -a'H^^^^ - [U + 2L {a'^ + a^') Xj) v^^ 

+ (c + {2L^a'^5ij - 2Lh'x, - AL"^ {a'^ + a^') x,Xj)) v. 

When ip is as in ()A.6p , one can see that A will not obey the hypotheses of Theorem 12.201 in 
particular the condition (j2.14p . unless L = 0. □ 

Appendix B. A weighted Sobolev inequality due to V. Maz'ya 

The following Sobolev inequality and corollary provide an alternative insight into Corollary 
18.61 and why it should not be expected to hold if the weights x^~^ and x^~^~^'^ are replaced by 
more general weights. Let Ap, 1 < p < oo, denote the class of Muckenhoupt weights in the sense 
of HH Definition 1.2.1 & 1.2.2]. 

Theorem B.l. [TH Theorem 2.6.1] Let w G Ai and let 1 < q < oo. Suppose that fi is a positive 
Radon measure, satisfying 

M:= sup V ' / s <oo. (B.l) 

Then there is a constant C = CM , where C only depends on d, q and the Ai constant for w, 
such that the following inequality holds for all u G C^iW^), 

\W\\Li{m.'l,p.) < C\\Du\\l^n^^d,^J,y (B.2) 

Conversely, if there exists a constant C such that ()B.2p holds for every u G Cq^{W^), then 
C > Cm , with C as before and, in particular, M is finite. 

Corollary B.2. Assume the hypotheses of Theorem I -B.il hold with 1 < q < oo. If q > 1, 

1 < p < q/{q — 1), and q < s < oo is defined by s = pq/{p — q{p — 1)), then for all u G C^{W^), 

ll'"llL«{R<*,/i) < C'i||Du||^p(iRd^^), (B.3) 
where Ci = pC /{p — q{p — 1)) and C is as in Theorem \B.1\ 

Proof. When p = 1, then ()B.3P follows from ()B.2p with s = g, so we may suppose p > 1. For 
1 < < oo defined by l/p + l/p' = 1, so p' = p/{p — I), and 7 > 1, we have 

\\\u\'^\\li{R'',,,) < 7C|||^^r~^^l^^lllLi(R^M) 

< 7C|||^^r"^llLP'(Rd,;,)ll^^^llLP(R^/.)- 

Now choose 7 to obey 79 = s = (7 — l)p', so that 

7 = — and s = — . 

P ~ 1\P ~ Ij P ~ Q{P ~ 1) 
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The preceding inequality becomes 

and this yields (|R3]l . □ 

It is useful to check Theorem IB . 1 1 when w = 1 and fi = dx. Then fi{Br{a)) = w{Br{a)) ~ r'^, 
so M in ()B.2p is finite if and only if sup^.^QV^'^'^/^'^'^ < oo, so 1 — d{l — 1/q) = and thus 
q = d/{d — I), to give the usual Sobolev inequality (compare [iQl Theorem 7.10]) when p = 1 and 
hence when I < p < d and d/{d — 1) < s < oo defined by s = dp/ (d — p). 

Note that to have s > 2 in ()B.3p . we need p > 2q/{3q — 2). If we also impose a constraint 
p < 2, then we must have q > 1, which holds in the basic case w = 1 and = dx. 



Appendix C. Weak maximum principles of Fichera, Oleinik, and Radkevic and the 

ELLIPTIC HeSTON OPERATOR 

We compare the weak maximum principles and uniqueness theorems provided by our article 
with those of Fichera, Olemik, and Radkevic [63] in the case of the elliptic Heston operator, A, 
in Example 11.21 on a domain 0' and show that those of Fichera, Olemik, and Radkevic are 
strictly weaker when < /3 < 1. 

C.l. Verification that the Heston operator and bilinear map coefficients have the 
required properties. We shall first illustrate how to choose tn so that (|8.2p holds for the 
elliptic Heston operator. A, on ^M. Denoting {x,y) = {xi,X2), we have 

Au = [uxx + 2gauxy + cr'^Uyy) - {r - q - y/2)ux - k{9 - y)uy + ru 

= iiyux + ygcruy)^ + {ygaux + ya^Uy)^ 

+ ^ {Q<yux + o-'^Uy) - {r - q- y/2)ux - n{0 - y)uy + ru 

= [{yux + yQcrUy)x + {yQcrUx + ya^Uy)^ 

+ {Qa/2 -{r-q- y/2)) Ux + [a'^ /2 - k(0 - y)) Uy + ru 
= -a'^Ux^Xj - b^Uxi + cu. 

Hence, 

2\g(7 a^J' V -(T^ + k{0 - y) y' 

As in [12], we choose 

1?= |(1- |p|)min{l,(T2} and « = /-^e"^l^l-'^3', 

so that log ro = (/3 — 1) log y — 'y\x\ — fj.y and 

(log Xo)x = -7 sign(2;) and (log Xo)y = (/3 - l)y~^ - fi. 

The coefficient matrix, a, weight w, and degeneracy coefficient, i?, clearly obey the regularity 
conditions (frajl . (frajl . (frajl . and (fCT]) . 
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Recalling that /? = IkQ ja^ and /i = 2k/ u"^ and recalling the definition (j8.24p of the coefficients 
6*, we see that 



while 



= h^ - ??(logrD)^a" - i9(log tt))ya^2 

= -Qa/2 + {r-q- y/2) + -ysign(x) - -(/3 - 1 - ^iy)Qa 
= {r -q- y/2) + ^ysign(x) - — {6 - y) 

gK6\ f 1 ■ , \ Qi^ 1 
- g + y o sign(x) H 



I? =b^ - ??(log lt))^.a^^ - ??(logtu)ya^ 



-cjV2 + - y) + |ysign(x)^a - ^(/? - 1 - ^y)^^ 
-ysign(x)^cr. 



The resulting bilinear map agrees with that in |181 Definition 2.2]. By making use of an affine 
change of variables |18^ Lemma 2.2], we may assume that r — q — qk9 /a = and so the expression 
for 6^ simplifies to 



b =y^-sign(x) + ---J. 



We can easily see that the coefficients, (a, b, c, d), of the bilinear map, a, associated with the elliptic 
Heston operator now obey the conditions (|7.13p . (j8.2p . (j8.3p . and ()7.12p : note that (d^) = 0. 

C.2. Comparison with the weak maximum principles and uniqueness theorems of 
Fichera. In the framework of Fichera (see [631 p. 308]), we lelF1 S denote the subset of points 
X £ dff where a^^ {x)ninj = (with n denoting the mwartf-pointing unit normal vector field, as 
in [631 p. 308]) and the Fichera function [63^ Equations (1.1.2) & (1.1.3)] (taking into account 
our sign convention in (II. 3p for the coefficients (a, 6, c) of A) is 



{V" - og) Uk = {b^ + {\ogxo)^^a^i^ Hk. 



Following [Ml P- 308], we denote by Si C S the subset where 6o > 0, by E2 C S the subset 
where 60 < 0) and by Sq C S the subset where 60 = 0; the set dff \ S is denoted by S3. By [63| 
Theorem 1.1.1], the characterization of the subsets S, Sq, Si, S2, S3 of the boundary dff remains 
invariant under smooth changes of the independent coordinates, (xi, . . . , Xri). 
In our example, we have S = dff n dM and n = (0, 1) along S, so thalo 

6o(x,0) = b\x,0) + {logn)^^a^^ = b\x,0) 
= _^2/2 + K^ = ^(/3-l). 



^'^In the work of Fichera [351 1581 [63l 164) . the boundary of the domain i^T" C R'' is usually denoted by E and 
is the subset of points a; G E where a^-' {x)ninj = 0. 

l^From [031 P- 310], one has that dff n dM is given by y = and -Ay = + ^.{6 ~ y) = bo, and thus 

T,'i = Ei for i = 0, 1, 2 in the notation of [63 p. 310]. 
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Hence, 

'Sa ifO</3<l, 



Si if/?>l, 
So if/3 = l, 



while E3 = Hrif?^. Th.e, first boundary value problem of Fichera \63\ Equations (1.1.4) & (1.1.5)] 
for the operator A is to find a function u G C'^{ff) such that 

= f on 0', ti = 5^ on S2 U S3, 

given a source function / on ^ and a boundary data function (7 on S2 U S3. But 



S2US3 



if < /3 < 1, 
n (9i^ if /3 > 1. 



Thus, for the Heston operator, the first boundary value problem of Fichera becomes 

Au = / on u = g on 



ii0</3 <1, 
ndff if/3>l. 



Therefore, we see that the first boundary value problem of Fichera differs from the formulations in 
[IHl [IHl EQI [30] when < /3 < 1, where a Dirichlet boundary condition along dffCidll is replaced by 
the requirement that u have a regularity property, in a weighted Holder or Sobolev sense, up to the 
boundary portion dff dM which is strictly weaker than that of the Fichera maximum principles 
[63^ Theorem 1.1.2 and Theorems 1.5.1 & 1.5.5] for C^(^) or H^^^{i^) functions, respectively. 
Note that 



SoUSi 



ifO</3<l, 
dmnd^ if/3>l. 



In the case of C^(^) functions on bounded subdomains 0" C H, we note that the Fichera maximum 
principle for C'^{ff) functions [Ml Theorem 1.1.2] requires that u G C2((^U Sq U Si) n C(^) and 
Au = f on ^USqUSi, which is stronger than the hypothesis of our Theorem 15.11 when /3 > 1, 
and yields, for r > 0, 

||n||c(^) < -||/llc(^?) V |bllc(S2UE3), 

whereO S2 U S3 = M n 9^ when /3 > 1 and S2 U S3 = dff when < /3 < 1. We see that the 
uniqueness result, when / = 0on ^USqUSi, afforded by the Fichera maximum principle [63[ 
Theorem 1.1.2] is weaker than that of our Theorem 15.11 when < /? < 1, since we only require 
5 = on HI n dff, and not g = on 30 to ensure that n = on ^. Indeed, the prescription of a 
Dirichlet boundary condition along 9IHI n 9^, when < /3 < 1, ensures that solutions to the first 
boundary value problem of Fichera are at most continuous up to dM. n dG and not smooth as in 

[ig [13 [101 [30]. 

Similar remarks apply to the Fichera maximum principle for weak solutions in L°°(^) [631 
Theorem 1.5.1 & 1.5.5]. Furthermore, our notions of weak solution, subsolution, or supersolution 
differ from those of [631 p. 318], which uses the adjoint operator A* to define these concepts for 



There is a typographical error in the statement of |63l Theorem 1.1.2], where E2 n E3 should be replaced by 
E2 U E3; compare [58] Theorem 1.1.2] 
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functions u £ L°°(^), with a bilinear map {u, A*v)i2(^ff^ and space of test functions v £ C" 
with 



= on 



if < /? < 1, 
ndff if/3>l, 



and thus imphes a Dirichlet boundary condition along dM D dff, when < /3 < 1, which is 
redundant in our framework of weighted Holder or Sobolev spaces. 



References 

1] A. Abdelkoddousse and A. Aurelien, Exact and high order discretization schemes for Wishart processes and 
their affine extensions, (2010), arxiv.org/abs/1006.2281 

M. Abramovitz and I. A. Stegun, Handbook of mathematical functions, Dover, New York, 1972. 
R. A. Adams, Sobolev spaces. Academic Press, Orlando, FL, 1975. 

R. F. Anderson, Diffusions with second order boundary conditions. I, Indiana Univ. Math. J. 25 (1976), 367- 
395. 

, Diffusions with second order boundary conditions. II, Indiana Univ. Math. J. 25 (1976), 403-441. 

S. R. Athreya, M. T. Barlow, R. F. Bass, and E. A. Perkins, Degenerate stochastic differential eguations and 
super-Markov chains, Probab. Theory Related Fields 123 (2002), 484-520. 

R. F. Bass and A. Lavrentiev, The submartmgale problem for a class of degenerate elliptic operators, Probab. 
Theory Related Fields 139 (2007), 415-449. 

R. F. Bass and E. A. Perkins, Degenerate stochastic differential equations with Holder continuous coefficients 
and super-Markov chains. Trans. Amer. Math. Soc. 355 (2003), 373-405. 

A. Bensoussan and J. L. Lions, Applications of variational inequalities m stochastic control, North-Holland, 
New York, 1982. 

S. Bonalede, A weak maximum principle and estimates o/esssupQ u for nonlinear degenerate elliptic equations, 
Czechoslovak Math. J. 46 (1996), 259-269. 

M. Borsuk, Second-order degenerate elliptic boundary value problems in nonsmooth domains, Sovrem. Mat. 
Fundam. Napravl. 13 (2005), 3-137. 

M.-F. Bru, Wishart processes, J. Theoret. Probab. 4 (1991), 725-751. 

J. Cox, J. IngersoU, and S. Ross, An intertemporal general equilibrium model of asset prices, Econometrica 53 
(1985), 363-384. 

, A theory of the term structure of interest rates, Econometrica 53 (1985), 385-407. 

M. G. Crandall, Viscosity solutions: a primer. Viscosity solutions and applications (Montecatini Terme, 1995), 
Lecture Notes in Math., vol. 1660, Springer, Berhn, 1997, pp. 1-43. 

M. G. Crandall, H. Ishii, and P-L. Lions, User's guide to viscosity solutions of second order partial differential 
equations. Bull. Amer. Math. Soc. (N.S.) 27 (1992), 1-67. 

C. Cuchiero, D. Filipovic, E. Mayerhofer, and J. Teichmann, Affine processes on positive semidefinite matrices, 
Ann. Appl. Probab. 21 (2011), 397-463. 

P. Daskalopoulos and P. M. N. Feehan, Existence, uniqueness, and global regularity for variational in- 
equalities and obstacle problems for degenerate elliptic partial differential operators in mathematical finance, 
arxiv.org/abs/1109. 1075vl 

P. Daskalopoulos and R. Hamilton, -regularity of the free boundary for the porous medium equation, J. 
Amer. Math. Soc. 11 (1998), 899-965. 

P. Daskalopoulos and E. Rhee, Free-boundary regularity for generalized porous medium equations, Commun. 
Pure Appl. Anal. 2 (2003), 481-494. 

J. Dolbeault, M. J. Esteban, M. Loss, and G. Tarantello, On the symmetry of extremals for the Caffarelli- 
Kohn-Nirenberg inequalities, Adv. Nonlinear Stud. 9 (2009), no. 4, 713-726. 

D. Dufiie, Credit risk modeling with affine processes. Journal ol Banking and Finance 29 (2005), 2751-2802. 
D. Dufiie, D. Filipovic, and W. Schachermayer, Affine processes and applications m finance, Ann. Appl. 
Probab. 13 (2003), 984-1053. 

D. Dufiie, J. Pan, and K. Singleton, Transform analysis and asset pricing for affine jump diffusions, Econo- 
metrica 68 (2000), 1343-1376. 

D. J. Duffy, Finite difference methods in financial engineering, Wiley, New York, 2006. 



PARTIAL DIFFERENTIAL OPERATORS WITH NON-NEGATIVE CHARACTERISTIC FORM 



65 



E. Ekstrom and J. Tysk, Boundary conditions for the single-factor term structure equation, Ann. Appl. Probab. 
21 (2011), 332-350. 

C. L. Epstein and R. Mazzeo, Degenerate diffusion operators arising in population biology, (2011), 341 pages, 
|arxiv . o rg/abs/1110 . 0032 

L. C. Evans, Partial differential equations, American Mathematical Society, Providence, RI, 1998. 

P. M. N. Feehan, Partial differential operators with non-negative characteristic form on stratified manifolds, 

maximum principles, and uniqueness theorems, Internet, math.rutgers.edu/~feehan 

P. M. N. Feehan and C. A. Pop, Degenerate-parabolic partial differential equations with unbounded coefficients, 
martingale problems, and a mimicking theorem for ltd processes, |arxiv.org/abs/1112.4824vl 

, Holder continuity for degenerate elliptic variational equations and inequalities in mathematical finance, 

I arxiv.o rg/abs/1110.5594vl 

, Stochastic representation of solutions to stationary and evolutionary degenerate variational equations 

and inequalities with Dirichlet boundary conditions, |arxiv.org/abs/1204. 1317vl[ 

W. Feller, Two singular diffusion problems. Annals of Mathematics 54 (1951), 173-182. 

G. Fichera, Sulle equazioni differenziali lineari ellittico-paraboliche del secondo ordine, Atti Accad. Naz. Lincei. 
Mem. CI. Sci. Fis. Mat. Nat. Sez. I. (8) 5 (1956), 1-30. 

, On a unified theory of boundary value problems for elliptic-parabolic equations of second order. Bound- 
ary problems in differential equations, Univ. of Wisconsin Press, Madison, 1960, pp. 97-120. 

D. Filipovic, Time-inhomogeneous affine processes. Stochastic Process. Appl. 115 (2005), 639-659. 

D. Filipovic and E. Mayerhofer, Affine diffusion processes: theory and applications. Advanced financial mod- 
elling. Radon Ser. Comput. Appl. Math., vol. 8, Walter de Gruyter, Berhn, 2009, pp. 125-164. 
J. Da Fonseca, M. Grasselli, and C. Tebaldi, A multifactor volatility Heston model. Quant. Finance 8 (2008), 
591-604. 

A. Friedman, Variational principles and free boundary problems, Wiley, New York, 1982, reprinted by Dover, 
New York, 2010. 

D. Gilbarg and N. Trudinger, Elliptic partial differential equations of second order, second ed., Springer, New 
York, 1983. 

A. Gnoatto and M. Grasselli, The explicit Laplace transform for the Wishart process, arXiv.org/abs/1107.2748. 
M. Grasselli and C. Tebaldi, Stochastic Jacobian and Riccati ODE m affine term structure models, Decis. 
Econ. Finance 30 (2007), 95-108. 

, Solvable affine term structure models. Math. Finance 18 (2008), 135-153. 

Q. Han and F. Lin, Elliptic partial differential equations, second ed., Courant Lecture Notes in Mathematics, 
vol. 1, Com'ant Institute of Mathematical Sciences, New York, 2011. 

S. Heston, A closed-form solution for options with stochastic volatility with applications to bond and currency 
options. Review of Financial Studies 6 (1993), 327-343. 

D. Kinderlehrer and G. Stampacchia, An introduction to variational inequalities and their applications. Aca- 
demic, New York, 1980. 

H. Koch, Non-Euclidean singular integrals and the porous medium equation, Habilitation Thesis, University of 
Heidelberg, 1999, www.mathematik.uni-dortmund.de/lsi/koch/publications.html 

N. V. Krylov, Lectures on elliptic and parabolic equations in Holder spaces, American Mathematical Society, 
Providence, RI, 1996. 

A. Kufner and B. Opic, Hardy type inequalities. Pitman Research Notes in Mathematics, Wiley, London, 1990. 
P. Lancaster and M. Tismenetsky, The theory of matrices, second ed.. Academic Press Inc., Orlando, FL, 1985. 
Ki-Ahm Lee, The obstacle problem for Monge-Ampere equation. Comm. Partial Differential Equations 26 
(2001), 33-42. 

G. M. Lieberman, Pointwise estimates for oblique derivative problems in nonsmooth domains, J. Differential 
Equations 173 (2001), 178-211. 

K. Miller, Barriers on cones for uniformly elliptic operators. Aim. Mat. Pura Appl. (4) 76 (1967), 93-105. 
D. D. Monticelli and K. R. Payne, Maximum principles for weak solutions of degenerate elliptic equations with 
a uniformly elliptic direction, J. Differential Equations 247 (2009), 1993-2026. 
J. R. Munkres, Topology: a first course, Prentice-Hall, Englewood Cliffs, N.J., 1975. 

N. Nadirashvili, A lemma on the inner derivative and the uniqueness of the second boundary value problem for 
second order elliptic equations, Dokl. Akad. Nauk. SSSR 261 (1981), 804-804, (Russian); English translation 
in Soviet Math. Dokl. 24 (1981), 598-601. 



66 



PAUL M. N. FEEHAN 



J. K. Oddson, On the boundary point principle for elliptic equations in the plane, Bull. Amer. Math. Soc. 74 
(1968), 666-670. 

O. A. Olemik and E. V. Radkevic, Second order equations with nonnegative characteristic form, Plenum Press, 
New York, 1973. 

M. A. Pozio, F. Punzo, and A. Tesei, Criteria for well-posedness of degenerate elliptic and parabolic problems, 
J. Math. Pures Appl. (9) 90 (2008), 353-386. 

P. Pucci and J. Serrin, The maximum principle. Progress in Nonlinear Differential Equations and their Appli- 
cations, 73, Birkhauser Verlag, Basel, 2007. 

F. Punzo and A. Tesei, On the refined maximum principle for degenerate elliptic and parabolic problems. 
Nonlinear Anal. 70 (2009), 3047-3055. 

, Uniqueness of solutions to degenerate elliptic problems with unbounded coefficients, Ann. Inst. H. 

Poincare Anal. Non Lineaire 26 (2009), 2001-2024. 

E. V. Radkevic, Equations with nonnegative characteristic form. I, J. Math. Sci. 158 (2009), 297-452. 

, Equations with nonnegative characteristic form. II, J. Math. Sci. 158 (2009), 453-604. 

J-F. Rodrigues, Obstacle problems in mathematical physics, North-Holland, New York, 1987. 
S. E. Shreve, Stochastic calculus for finance. Volume II: Continuous-time models. Springer, New York, 2004. 
D. W. Stroock and S. R. S. Varadhan, Multidimensional diffusion processes, Springer, Berlin, 1979. 
K. Taira, Semigroups, boundary value problems and Markov processes, Springer, Berlin, 2004. 

G. M. Troianiello, Elliptic differential equations and obstacle problems. Plenum Press, New York, 1987. 
N. S. Trudinger, Maximum principles for linear, non-uniformly elliptic operators with measurable coefficients, 
Math. Z. 156 (1977), 291-301. 

B. O. Turesson, Nonlinear potential theory and weighted Sobolev spaces. Lecture Notes in Mathematics, vol. 
1736, Springer, Berlin, 2000. 

A. Zettl, Sturm Liouville theory, American Mathematical Society, Providence, RI, 2005. 

R. Zvan, P. A. Forsyth, and K. R. Vetzal, Penalty method for American options with stochastic volatility, 
Journal of Computational and Applied Mathematics 91 (1998), 199-218. 

Department of Mathematics, Rutgers, The State University of New Jersey, 110 Frelinghuysen 
Road, Piscataway, NJ 08854-8019 

E-mail address: f eehanOmath. rutgers . edu 



